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SUMMARY 


Tropospheric  bias  errors  are  highly  predictable  using 
radiosonde  or  refractometer  profiles;  residual  errors  from 
1%  to  3%  of  the  initial  bias  levels  are  commonly  attained 
using  procedures  described  in  Section  3 .  Data  to  within  one- 
half  foot  in  range  and  20  to  70  //radians  in  angle  can  be  ex¬ 
pected  at  elevation  angles  above  five  degrees . 

Tropospheric  fluctuation  errors  are  not  correctable 
using  any  known  procedure,  and  will  cunount  to  a  few  tenths 
of  a  foot  in  range,  and  10  to  50  //radians  in  angle  (depend¬ 
ing  on  the  baseline  or  aperture  used  for  measurement) ,  under 
normal  weather  conditions . 

The  relationship  between  temporal  and  spatial  correlation 
of  tropospheric  fluctuations  has  been  investigated,  based  on 
data  obtained  by  the  National  Bureau  of  Standards.  The  effect 
of  short -period  fluctuations  is  described  by  Figs.  2-7  and 
2-8,  and  is  consistent  with  a  drift  of  tropospheric  anomalies 
across  a  fixed  measurement  path  at  the  speed  of  the  prevailing 
wind. 

In  range  instrumentation  applications,  where  the  beam  is 
not  fixed,  the  residual  "bias"  and  long-term  error  components 
will  change  as  the  beam  moves,  and  additional  atmospheric 
rate  errors  will  be  generated,  as  shown  in  Fig.  4-6.  These 
errors  will  be  proportional  to  the  tangential  velocity  of 
the  missile,  and  will  typically  be  five  to  fifty  times  the 
errors  measured  for  a  fixed  beam. 

The  uncertainty  in  tropospheric  path  leads  to  errors 
equivalent  to  motion  of  the  instrument  on  the  ground.  The 
motion  of  the  "virtual  source"  typically  amounts  to  several 
feet  normal  to  the  path  and  a  few  tenths  of  a  foot  along  the 
path . 
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Ionospheric  errors  are  essentially  unpredictable,  and  will 
exceed  the  residual  tropospheric  errors  when  operating  fre¬ 
quencies  below  3000  mcps  are  used.  Even  in  the  5000-6000  mcps 
band  the  ionospheric  errors  will  contribute  to  overall  atmos¬ 
pheric  error  during  daytime  operation. 

Redundant  measurements  performed  at  two  frequencies  below 
3000  mcps  can  be  used  to  correct  for  ionospheric  error  in 
both  range  and  angle. 

The  lowest  atmospheric  errors  are  found  in  trilateration 
systems  using  very  long  baselines.  Total  position  and  velocity 
errors  for  a  typical  satellite  track  (660  miles  range,  100 
miles  altitude)  through  average  weather,  are  as  follows: 


RMS  Position 

RMS  Velocity 

Error  (feet) 

Error  (ft/sec) 

Range -angle  tracker  at 

310 

16 

6000  mcps 

Interferometer , 

100 

2.4 

10 , 000  mcps  (Mistram) 

Trilateration  system. 

2000  mcps 

19 

0.9 

6000  mcps 

2.5 

0.15 

The  above  errors  may  be  increased  or  decreased  by  a 
factor  of  two  or  three  for  different  weather  conditions 
(and  a  2000  mcps  for  different  ionospheric  conditions) . 

The  trilateration  errors  shown  are  dependent  upon  perfect 
survey  of  station  location,  as  well  as  instrumental  errors 
below  one -ha If  foot  in  range  and  0.02  ft/sec  in  range -rate. 


1 .  INTRODUCTION 


The  Ad  Hoc  Panel  on  Electromagnetic  Propagation  was 
convened  as  part  of  a  continuing  effort  by  the  National  Academy 
of  Sciences -National  Research  Council  on  behalf  of  the  Atlantic 
Missile  Range  as  requested  by  Headquarters ,  Air  Force  Systems 
Command . 

The  Ad  Hoc  Panel  on  Basic  Measurements  discussed  in  its 
Report  how  well  we  can  measure  at  the  present  time,  such 
fundamental  quantities  as  a  length,  time,  and  the  velocity  of 
light .  The  Report  notes  that  there  are  fundamental  limita¬ 
tions  to  tracking  accuracy  imposed  by  our  inability  to  measure 
these  basic  quantities  with  more  precision,  but  it  points  out 
that  state  of  the  instrumentation  art  has  not  yet  approached 
these  limitations.  There  are,  however,  other  fundamental 
limitations  to  tracking  accuracy  which  today  are  beginning  to, 
if  not  already,  restricting  the  capability  of  our  instrumenta¬ 
tion  state-of-the-art.  One  of  these  major  limitations  is  that 
imposed  by  atmospheric  refraction.  Consequently,  the  Ad  Hoc 
Panel  on  Electromagnetic  Propagation  was  formed  to  consider 
this  problem. 

This  Panel  met  on  11  May  1962  in  Washington,  D.  C.  and 
heard  discussions  as  to  how  tracking  accuracy  requirements 
were  arrived  at  for  one  particular  program,  of  the  current 
tracking  capabilities  of  the  AMR,  and  of  various  research  work 
which  instrumentation  and  atmospheric  physics  people  are  con¬ 
ducting.  Due  to  the  quantity  and  divergence  of  the  material 
presented,  the  Panel  could  not  arrive  at  a  consensus  of  opinion 
regarding  a  report.  Consequently,  Mr.  David  K.  Barton  was 
appointed  Editor  of  the  Panel's  Report  by  the  chairman.  Mr. 
Barton  drafted  the  Report  of  the  Panel  from  material  contri¬ 
buted  by  Dr.  Robert  S.  Fraser,  Dr.  John  B.  Smyth,  Mr.  Preston 
Landry,  and  himself. 

This  Report  of  the  Ad  Hoc  Panel  on  Electromagnetic  Propa¬ 
gation  has  attempted  to  state  the  current  extent  of  our  know¬ 
ledge  concerning  atmospheric  refraction  and  its  effect  on 
tracking  accuracy,  which  should  be  of  as  much  benefit  to  the 
range  users  as  to  the  range  operators ,  and  also  has  made 
recommendations  which,  if  followed,  should  lead  to  increased 
accuracy  both  in  the  near  and  in  the  more  distant  future. 
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2.  CLASSIFICATION  AND  DESCRIPTION  OF  ATMOSPHERIC  ERRORS 
2 . 1  Classification  of  Atmospheric  Errors  by  Type. 

There  are  many  ways  of  describing  propagation  errors  in 
precision  tracking  systems.  The  four  classifications  shown 
in  Table  I  are  suggested  to  cover  the  characteristics  of 
most  interest  to  the  developers  and  users  of  missile  and 
space  range  instrvimentation : 

Table  I 


a.  Source  of  Error: 

Tropospheric 

Ionospheric 

b.  Measured  Quantity 

Angle  of  arrival  or  phase  difference 
Range  delay  or  signal  phase 

c.  Spatial  Correlation  of  Error: 

Across  radar  aperture  (5  to  100  feet) 

Across  short  baseline  (100  to  1000  feet) 
Across  long  baseline  (1000  to  100,000  feet) 

d.  Temporal  Correlation  of  Error: 

Bias  (fixed  during  one  track) 

Fluctuation  (periods  of  seconds  or  minutes) 


For  each  combination  of  the  above  characteristics,  the 
error  should  be  known  as  a  function  of  operating  frequency, 
target  altitude,  elevation  angle  (or  slant  range)  and  state 
of  the  atmosphere.  Except  in  rare  instances,  the  instru¬ 
mentation  system  may  be  assumed  to  be  at  sea  level. 

2.2  Effects  of  the  Lower  Atmosphere 

The  propagation  errors  caused  by  the  lower  atmosphere 
will  be  analyzed  in  this  section.  The  lower  atmosphere  for 
the  purposes  of  this  discussion  will  be  defined  as  the  atmos¬ 
phere  below  a  height  of  40  km.  The  object  of  the  analysis 
is  to  find  the  accuracies  that  the  position  and  motion  of  a 
target  above  the  lower  atmosphere  can  be  measured  by  means 
of  radio  systems  on  the  surface  of  the  earth.  The  propaga¬ 
tion  effects  can  not  be  presented  completely  in  this  brief 
discussion.  Hence,  the  location  and  motion  of  the  target 
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will  be  restricted,  in  order  to  demonstrate  simply  the  propa¬ 
gation  effects  encountered  most  frequently. 

a.  Refractive  Index  of  Air 

The  refractive  index  of  any  substance  is  composed  of 
a  real  and  an  imaginary  part.  The  phase  velocity  of  an  electro¬ 
magnetic  wave  depends  on  the  real  part.  The  attenuation  of  the 
electromagnetic  energy  depends  on  the  imaginary  part.  The 
attenuation  is  important  to  radio  tracking  and  guidance  systems , 
if  the  radio  signal  becomes  too  weak  for  the  system  to  operate 
on  it.  The  attenuation  by  the  atmospheric  gas  is  negligible 
for  most  purposes  at  frequencies  slightly  below  22,000  mc/s, 
which  is  the  lowest  frequency  of  an  absorption  band  of  the 
atmospheric  gas.  Atmospheric  particulate  matter  such  as  rain, 
snow,  and  clouds  also  attenuate  electromagnetic  energy  passing 
through  them.  These  latter  effects  are  rarely  important  for 
many  tracking  and  guidemce  systems.  Consequently,  this  dis¬ 
cussion  of  propagation  errors  will  apply  only  to  radio  fre¬ 
quencies  less  than  20,000  mc/s,  where  atmospheric  attenuation 
is  not  important  for  many  tracking  and  guidance  systems . 

The  real  part  of  the  index  of  refraction  (n)  of  the 
atmospheric  gas  can  be  computed  from  the  expression 

(n  -  1)  .  10®  =  N  =  (P  +  iSlOjj  )  ^  (2_i) 

T  T 

where  T  represents  the  temperature  of  the  atmosphere  in  degrees 
Kelvin;  P  represents  the  total  pressure  of  the  atmosphere  in 
millibars;  p  represents  the  partial  water  vapor  pressure  in 
millibars;  and  N,  which  is  called  the  refractivity ,  is  used 
frequently  in  radio  propagation  work  instead  of  the  index  of 
refraction.  The  refractivity  can  be  calculated  from  Eq. (2-1) 
with  an  accuracy  of  0.5  per  cent  for  the  usual  range  of  atmos¬ 
pheric  variables  and  for  radio  frequencies  less  than  30,000 
mc/s.  Since  only  radio  frequencies  less  than  20,000  mc/s  will 
be  considered  in  this  report,  the  propagation  effects  to  be 
discussed  will  be  independent  of  radio  frequency. 

The  dependence  of  the  refractivity  N  on  height  is 
shown  in  Fig. 2-1.  Mean  vertical  profiles  for  summer,  winter, 
and  an  elevated  station  (Ely)  are  shown.  Also,  the  approximate 
extreme  values  are  indicated.  The  values  of  N  near  sea  level 
vary  from  about  270  to  405.  The  mean  value  of  N  near  sea  level 
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Altitude  Atxxte  Mean  Sea  Level  in  Kilometere 


Figure  2-1.  Means  and  extremes  of  radio-re fractlvlty 
height  structure,  and  a  comparison  with  extreme 
4/3  earth  profiles  (Bean  and  Thayer,  1959). 
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is  313 .  The  comparatively  large  difference  in  N  decrease 
with  increasing  height .  The  value  of  N  is  about  one  at  a 
height  of  40  km.  Later,  use  will  be  made  of  the  fact  that 
the  N-vertical  profile  can  be  approximated  by  an  exponential 
function  with  sufficient  accuracy  to  estimate  many  propagation 
errors . 

b.  Range  and  Angle  Errors  for  a  Mean  Atmosphere 

Because  the  index  of  refraction  of  the  earth's  lower 
atmosphere  is  greater  than  one,  the  apparent  position  of  a 
target  which  is  determined  by  a  radio  system  differs  from  the 
position  which  the  S2une  system  would  measure  if  there  were  no 
atmosphere.  The  apparent  position  measured  in  the  presence 
of  an  atmosphere  minus  the  true  position  is  the  error  in  the 
position  measurement.  Important  propagation  effects  are 
revealed  when  the  position  errors  are  calculated  for  a  simple 
model  of  the  earth's  atmosphere,  such  as  one  that  is  spherically 
stratified.  In  other  words,  the  refractivity  N  for  such  a 
model  is  constant  on  an  arbitrary  spherical  surface  which  is 
concentric  with  the  center  of  a  spherical  earth.  Furthermore, 

N  does  not  change  with  respect  to  time. 

Range  Errors 

The  range  error  measured  by  a  radar  is  given  by  the 
expression 

4R,  =  /  n(z)  CSC  ^0(z)j  dz  -  Rq,  (2-2) 

o 

where  the  integral  represents  the  distance  measured  by  the 
radar  and  R^  represents  the  true  distance;  6(z)  is  the  eleva¬ 
tion  angle,  or  the  angle  between  the  curved  radio  ray  and  the 
spherical  surface  which  the  ray  intersects  at  height  z  above 
the  radar  (z  is  measured  normal  to  the  spherical  surfaces) ; 
and  h  is  the  distance  between  the  two  spherical  surfaces  that 
pass  through  the  radar  and  the  target.  The  integration  above 
a  height  of  70  km.  is  negligible  for  most  guidance  and  track¬ 
ing  purposes  since  the  index  of  refraction  is  very  small. 
(Ionospheric  effects  are  not  being  considered) .  If  the 
initial  elevation  angle  at  the  radar  is  large  enough.  Eg. (2-1) 
can  be  written  as 

6 

4Rg  *  CSC  Oq  /  N(z)  .  10  dz  +  /  CSC 

o  o 
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[e(z)l  dz  -  Rq  • 
(2-3) 


The  last  two  terms  on  the  right  of  the  above  equation  repre¬ 
sents  the  difference  between  the  length  of  the  curved  radio 
ray  and  the  true  distance.  If  the  initial  elevation  angle  is 
greater  than  10°  (0q  >  10°) ,  these  two  terms  may  be  neglected 
in  comparison  with  the  first  term  on  the  right  of  Eq.  (2-3) . 
The  latter  term  represents  the  apparent  increase  in  distance 
to  the  target,  because  the  phase  velocity  of  the  radio  wave 
passing  through  the  atmosphere  is  less  than  it  is  in  free 
space.  In  order  to  estimate  the  magnitude  of  this  term, 
assume  that  the  refractivity  N  decreases  exponentially  with 
increasing  height;  that  is. 


N{z) 


(2-4) 


where  Nq  is  the  value  of  the  refractivity  at  the  radar ;  H  is 
a  scale  height ;  and  z  is  the  height  above  the  radar .  When  this 
value  of  N(z)  is  substituted  into  Eq.  (2-3),  the  expression 
for  the  range  error  becomes 

i  N^H  CSC  .  10-6  (1  -  e-^/”)  .  (2-5) 


The  range  error  increases  as  the  path  length  through  the 
atmosphere  increases,  or  as  the  initial  elevation  angle  (0©) 
decreases.  To  compute  the  range  error  for  a  particular  case, 
use  the  following  values  of  the  parauneters :  the  average  value 
of  the  refractivity  at  sea  level,  N©  =  313;  the  average  scale 
height,  H  =  7  km;  the  height  of^the  target  h  =  30  km.;  and  the 
initial  elevation  angle  0^  =  20  .  The  range  error  for  this 
case  is  =  6.4m.  Detailed  computations  of  range  errors 
for  an  exponential  atmosphere  are  given  by  Bean  and  Thayer 
(1959) .  The  range  error  ^R©  computed  by  Bean  and  Thayer  for 
the  CRPL  Exponential  Reference  Atmosphere  with  Ng  =  313  is 
plotted  in  Pig.  2-2,  as  a  function  of  slant  range  R©  and 
elevation  angle  0©. 

Angular  Errors 

Next,  the  angular  errors  will  be  discussed.  The 
errors  which  appear  in  interferometer  measurements  of  angle 
will  be  derived  before  presenting  the  angular  errors  associ¬ 
ated  with  a  tracker.  The  following  derivation  of  the  inter¬ 
ferometer  errors  is  based  on  a  report  by  Thayer  and  Bean 
(1962) .  The  geometry  for  an  inter forometer  system  is  shown 
in  Figure  2-3.  Ai  and  A2  indicate  the  two  antennas  which  are 
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Figure  2-2 

Slant  Range  R^  in  Nautical  Miles 


Range  Bias  vs.  Range  for  CRPL  Exponential 
Reference  Atmosphere  (N^  *  313) 


Figure  2-2 


Figure  2-3.  Qeoaetry  of  an  Interferometer 


radio  system. 
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separated  by  a  distance  B.  Ri  and  R2  indicate  the  true 
distances  between  the  target  and  the  corresponding  antennas. 
Note  that  j3  '  refers  to  the  angle:  between  a  radio  ray  from 
antenna  1  to  the  target  and  the  baseline,  and  not  the  angle 
between  the  radio  ray  and  the  tangent  plane.  An  exact 

B  2R  4s  ’ 


solution  for  0 '  is 


=  sin 


-1 


1  - 


where 

dS  =  R^  -  R2  . 


In  order  to  show  simply  the  refraction  errors  which  depend 
on  a  spherically  stratified  atmosphere,  the  following 
assumptions  are  made : 


a. )  the  range  R  to  the  target  is  large 

compared  to^the  baseline  length  B, 
say  >100  B;  and 

o 

b. )  the  value  of  is  larger  than  10 


Then  an  expression  for  the  error  in  the  angle 
caused  by  a  spherically  stratified  atmosphere  is 


40' 


- 

B  sin 


d|3') 


(2-7) 


The  Seune  expression  for  the  error  A  /3'  results  from  the 
derivation  for  a  plane-parallel  atmosphere.  The  numerator 
on  the  right  of  Eq.  (2-7)  represents  the  difference  in  range 
errors  caused  by  the  atmosphere;  and  the  denominator  on  the 
right  of  Eq.  (2-7)  represents  the  effective  baseline  length, 
or  the  perpendicular  distance  between  two  rays  from  the 
distant  target  to  antennas  A  and  A  .  If  the  expression  for 
the  range  error  (2-3)  is  substituted  into  (2-7) ,  then  the 
expression  for  the  angular  error  can  be  approximated  by 

Ab'  =  f  N(z)  .  lo’^  dz.  >10°,  (2-8) 

•'o 
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where  is  the  angle  that  a  radio  ray  at  A  makes  with  the 
baseline  ?  3*  is  greater  than  10° ;  r^  is  the'^radius  of  the 
earth;  and  both  antennas  are  at  the  same  height.  The  refrac¬ 
tion  error  for  an  interferometer  is  not  zero,  even  when  the 
atmosphere  is  spherically  stratified.  The  error  is  nearly 
independent  of  the  distance  between  the  two  antennas  for  the 
assumptions  that  have  been  made.  The  error  is  proportional 
to  the  cube  of  the  cosec  ;  and  the  error  depends  on  the 
vertical  profile  of  the  refractivity  N.  However,  if  N  is 
assumed  to  decrease  exponentially  with  increased  height, 
then  (2-4)  can  be  substituted  into  (2-8)  to  obtain 


4  0'* 


3  B* 
csc^  ^ 

r 

o 


NoH- 


(1  -  e-^/») 


(2-9) 


If  the  target  is  sufficiently  high,  the  exponential  is  negli¬ 
gible,  and  the  refraction  error  is  directly  proportional  to 
the  refractivity  at  the  antennas. 

The  elevation  angle  error  for  a  tracker  has  a 
different  dependence  on  the  initial  elevation  angle.  If  the 
refractivity  (N)  decreases  exponentially  with  increasing 
height  and  if  the  Initial  elevation  angle  is  not  too  small, 
then  the  tracker  elevation  angle  error  is  given  by  the  ex¬ 
pression 

«  =  Nq  •  10“®  ctn  ©o  ,  (2-10) 


where  Bq  is  the  value  of  the  elevation  angle  at  the  tracker. 

The  tracker  elevation  angle  error  is  also  directly  propor¬ 
tional  to  the  refractivity  at  the  tracker  but  decreases  less 
rapidly  with  respect  to  decreasing  elevation  angle  than  does 
the  interferometer  elevation  angle  error.  The  angle  error 
computed  by  Bean  and  Thayer  for  the  CRPL  Exponential  Reference 
Atmosphere  with  =  313  is  plotted  in  Figure  2-4  as  a  function 
of  slant  range  and  elevation  angle  0q. 

Both  the  tracker  and  interferometer  elevation  angle 
refraction  errors  are  shown  for  mean  conditions  in  Figure 
2-5.  Both  errors  decrease  with  increasing  elevation  angle, 
but  the  interferometer  errors  decrease  more  rapidly  than  the 
tracker  errors,  except  for  large  elevation  angles.  The 
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Slant  Range  R^  In  Nautical  Milas 
10  30  100  300  1000 


30  100  300  1000  3000 

Figure  2-4  Slant  Range  R^  in  Kilometers 

Tracker  Elevation  Angle  Error  vs.  Range 
for  CRPL  Exponential  Reference  Atmosphere 


MEAN  ELEVATION  ANGLE  ERROR  IN  MICRORADIANS 

rigux*  2-5.  Htan  tracKar  and  IntarfaroaMtar 

rafractlon  arrora.  Tha  curva  for  8*3  kB  alao 
rapraaanta  tha  xaaidual  traekar  arror  aftar 
oorzactlon  according  to  Sq.  2-10. 


CYCLES  PER  SECOND 

Figura  2-6.  spactra  of  rafractlvity  and  apparant 
ranga  varlationa  at  Maul  and  Colorado  (Thonpaon, 
Janaa  and  Kirkpatrick  (1960). 
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tracker  error  is  greater  than  the  interferometer  error, 
except  at  small  elevation  angles  where  the  effective  base¬ 
line  length  of  the  interferometer  becomes  small.  The  inter¬ 
ferometer  errors  are  nearly  Independent  of  the  distance 
between  the  antennas  (B)  for  large  elevation  angles,  but  the 
errors  are  much  smaller  for  the  longer  baselines  at  the  small 
elevation  angles. 

The  interferometer  error,  for  the  shorter  baselines, 
is  approximately  equivalent  to  the  tracker  error  corrected 
by  Bq.  (2-10)  for  surface  refractivity.  In  such  an  inter¬ 
ferometer,  the  term  NQlO“°ctn  6©  is  supplied  by  the  added 
delay  in  the  path  to  antenna  after  the  wavefront  has 
reached  the  antenna  A..  Any  deviation  from  the  conditions 
that  apply  in  Eq.  (2-10)  will  cause  the  interferometer  data 
and  the  corrected  tracker  data  to  depart  from  the  curves 
shown  in  Figure  2-5. 

The  interferometer  azimuthal  refraction  errors  are 
negligibly  small  for  the  condition  that  the  two  rays  from 
the  target  to  the  two  zmtennas  are  nearly  parallel.  In  this 
case,  both  rays  are  refracted  the  same.  The  tracker  azimuthal 
error  Is  zero,  since  the  refraction  is  proportional  to  the 
gradient  of  the  index  of  refraction;  the  gradient  is  zero  in 
the  azimuthal  direction  for  a  spherically  stratified  atmosphere. 

Rate  Errors 

The  propagation  rate  bias  errors  of  a  guidance  system 
operating  in  a  spherically  stratified  atmosphere  can  be  easily 
derived  from  the  formulas  for  the  position  errors  by  taking 
derivatives  with  respect  to  time.  The  expression  for  the  range 
rate  error,  as  obtained  from  (2-5) ,  is 

i  -  ctn  Bq  4Rg,  h  >32  km.  (2-11) 

The  range  rate  errors  are  proportional  to  the  rate  of  change 
of  the  elevation  angle.  For  an  example,  let  6©  =  20° »  ®o  =10“ 
rad/sec,  amd  4  R^  =  6.4  m;  then  the  range  rate  error  is 
4  ftg  =  -  0.02  m/sec.  The  angular  rate  error  for  an  interfer¬ 
ometer  is  obtained  by  differentiating  (2-9)  with  respect  to 
time : 

*  h 

4*3'  =  -  3  ctn  ^  |3  41?.  k  >32  km.  (12-12) 
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As  an  example,  let  ^  be  the  elevation  angle  and  equal  to 
20°;  let  p  =  10“^  rad/sec;  then  =  0.06  |jir/sec. 

c.  Geographical,  Seasonal,  and  Diurnal  Variations 
in  Errors _ _ 

Corrections  can  be  made  for  the  propagation  errors 
that  are  caused  by  a  spherically  stratified  atmosphere.  As 
has  been  shown,  the  errors  are  nearly  proportional  to  the 
surface  value  of  the  refractivity  (Nq) ,  if  the  target  is  far 
enough  away  and  if  the  elevation  angle  at  the  antennas  exceeds 
about  10°.  The  range  and  refraction  corrections  can  be  based 
on  various  averages  of  the  surface  refractivity,  such  as  the 
annual,  seasonal,  or  diurnal  means,  or  on  the  value  of  the 
refractivity  at  the  radio  antennas  a  short  time  before  the 
radio  system  is  to  be  used  for  guidance  or  tracking.  Then 
the  standard  deviation  of  the  propagation  errors  are  pro¬ 
portional  to  the  standard  deviation  of  the  surface  refrac¬ 
tivity.  A  good  compilation  of  climatic  refractive  data  has 
been  published  by  Bean,  Horn,  and  Ozanich  (1960) .  The  standard 
deviation  of  the  surface  refractivity  has  been  estimated 
approximately  from  that  data  and  put  in  Table  2-2.  The  stand¬ 
ard  deviations  for  the  range  and  elevation  angle  errors  that 
appear  in  Table  1  have  been  computed  from  formulae  (2-5)  and 
(2-9) ,  respectively,  where  the  propagation  errors  and  the 
refractivity  in  those  formulae  were  replaced  by  the  standard 
deviations  of  the  corresponding  quantities.  The  standard 
deviation  of  the  refractivity  which  is  measured  one  hour  before 

TABLE  2-2.  Standard  deviation  of  range  and  inter¬ 
ferometer  elevation  angle  errors .  The  elevation 
angle  is  assumed  to  be  20^  and  the  height  of  the 
target  exceeds  32  km.  The  temporal  standard  devi¬ 
ations  of  the  refractivity  apply  to  Mieuni. 
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the  radio  system  is  used  to  guide  or  track  is  based  on  the  fact 
that  the  refractivity  can  be  measured  with  an  accuracy  of 
4  N-units  and  will  change  3  N-units  during  one  hour.  The  data 
in  Table  2-2  show  that  as  the  period  or  geographical  extent 
for  which  the  error  applies  decreases,  the  error  decreases. 

Bean  and  his  co-workers  have  developed  empirical 
methods  for  making  refraction  and  range  corrections.  The 
methods  involved  one  to  three  parameters,  the  most  Important 
being  the  surface  refractivity.  These  empirical  corrections 
reduce  the  propagation  errors  caused  by  a  spherically  strati¬ 
fied  atmosphere  to  a  few  per  cent  of  the  uncorrected  errors. 

d.  Random  Errors 


The  lower  atmosphere  is  not  quite  spherically 
stratified,  and  its  refractivity  is  constantly  changing. 

Hence,  the  propagation  errors  are  actually  larger  than  have 
been  given  for  a  spherically  stratified  atmosphere.  In 
general,  the  apparent  position  of  a  target  fluctuates  about 
its  true  position.  The  fluctuating  position  and  rate  errors 
depend  on  how  the  propagation  errors  pass  through  several 
filters  that  are  either  inherent  to  or  designed  in  any  guidance 
or  tracking  system. 

Hence ,  it  is  useful  to  review  how  the  propagation 
errors  are  filtered  by  a  radio  system.  Such  filters  are 
discussed  in  an  article  by  Wheelon  (1959) .  Consider  a 
propagation  error,  such  as  range  error  for  excunple,  to  be 
a  stochastic,  stationary  function  of  time  and  to  have  the 
spectral  representation 

oe  . 

x(t)  =  f  ds{w)  ,  (2-13) 

-oo 

where  W  =  2iTf,  f  being  the  frequency  of  the  radio  wave ,  EdS  (<«))  =0 , 
E  dS((i>)  dS*(u>)=  s((i;)d(t).  E  denotes  the  expectation  operator 
and  s  (cj)  the  spectral  density  of  the  process.  When  x(t)  is 
the  input  to  a  linear,  time-invariant,  passive  filter,  which 
has  the  frequenty  response  function  1 (W)  ,  the  spectral  density 
of  the  output  is  1 1  ((«»)|  ^  s  (w)  .  The  absolute  square  of  the 
frequency  response  function  for  rate  errors  is  which 
implies  that  the  rate  errors  are  nearly  independent  of  the 
low  frequency  errors.  The  finite  size  of  the  receiving 
f  antenna  acts  as  a  filter  which  does  not  pass  errors  caused 
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by  atmospheric  irregularities  smaller  than  the  cross  section 
of  the  antenna  aperature.  Aperture  filtering  is  not  important 
for  present  guidance  and  tracking  systems,  since  the  propa¬ 
gation  errors  at  frequencies  above  the  aperature  high  fre¬ 
quency  cutoff  are  generally  negligible  with  respect  to  the 
low  frequency  errors  which  are  passed.  Linear  filtering  of 
the  data  for  a  time  r  has  the  absolute  square  frequency 
response  function  of 


(sinioZi)  2  ,  (2_14) 

0)t/2 

This  filter  passes  low  frequency  errors.  Since  the  linear 
filter  has  a  high-frequency  cutoff  at  about  the  inverse  of 
the  smoothing  time  r  for  some  of  the  propagation  errors,  the 
vari2uice  of  these  errors  decreases  as  the  smoothing  time  in¬ 
creases  .  Some  radio  systems  use  much  more  sophisticated 
time  filters  than  linear  filters.  Another  filtering  action 
depends  on  the  procedure  for  applying  corrections  for  propa¬ 
gation  errors.  For  exaunple,  a  tracker  elevation  angle  re¬ 
fraction  correction  can  be  easily  computed  from  (2-10) .  If 
T  is  the  interval  of  time  between  the  measurement  of  the 
refractivity  which  is  used  in  the  correction  formula  and 
the  use  of  the  correction  by  the  radio  system,  then  the 
absolute  square  of  the  frequency  response  function  is 

4sin^(wT/2).  (2-15) 

If  T  =  0,  the  correction  is  not  in  error  (assigning  that  the 
correction  cem  be  determined  accurately) .  As  T  becomes  large, 
the  variance  of  the  error  approaches  twice  the  variance  which 
would  result  if  no  correction  were  applied.  This  filtering 
action  has  been  indicated  in  section  c  about  temporal  and 
geographical  variations  of  the  refractivity.  Besides  the 
corrections  for  the  low  frequency  propagation  errors,  which 
are  sometimes  considered  as  a  bias  correction,  the  filtering 
of  the  propagation  errors  in  time  is  the  most  important 
filtering  action  for  many  systems.  Hence,  the  effect  of  this 
filter  (2-14)  on  measured  errors  will  be  discussed  in  the 
following  discussion. 

1.  Maui  Experiments 

Valuable  data  on  radio  propagation  errors  have  been 
collected  from  experiments  made  at  Maui,  Hawaii,  by  the 
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National  Bureau  of  Standards.  The  experiments  were  designed 
to  measure  range  and  angle  errors  that  occurred  during  periods 
of  the  order  of  minutes.  The  propagation  path  extended  a 
distance  of  2 Sion,  from  near  sea  level  to  a  summit  3  Ion.  high. 
The  elevation  angle  of  the  propagation  paths  was  about  7°. 

The  interferometer  measvired  azimuthal  errors,  since  the  inter¬ 
ferometer  baseline  was  approximately  level  and  perpendicular 
to  a  line  from  one  antenna  to  the  transmitter  at  the  summit. 
The  propagation  path  extended  through  only  part  of  the  lower 
atmosphere,  but  that  path  was  chosen  because  it  was  believed 
that  large  propagation  errors  would  occur  there.  The  propa¬ 
gation  errors  observed  at  Maui  have  been  consistent  with 
propagation  errors  which  have  been  measured  while  tracking 
moving  targets  above  the  lower  atmosphere  and  at  elevation 
angles  greater  than  15°.  The  speed  at  which  such  targets 
were  moving  is  not  given  for  security  reasons.  A  svimmary  of 
the  propagation  errors  measured  at  Maui  follows.  The  data 
are  taken  from  a  Space  Technology  Laboratory  report  (1958) 
and  Norton  et  al  (1961) . 

Rauige 

Power  spectra  of  apparent  range  and  refractivity 
fluctuations  are  shown  in  Figure  2-6.  The  spectral  values 
for  frequencies  higher  them  10~  sec~  are  the  median  values 
of  many  spectra.  The  parts  per  million  indicated  on  the 
ordinate  refers  to  either  the  change  in  path  length  from  some 
arbitrary  reference  divided  by  the  path  length,  or  to  the 
refractivity  times  10”  .  The  Colorado  data  shown  apply  to  a 
horizontal  path  near  Boulder,  15  km.  long,  and  about  100  m. 
above  the  terrain.  The  spectral  density  of  range  errors  in¬ 
creases  with  decreasing  spectral  frequency  and  converges 
to  the  spectral  density  of  the  refractivity  fluctuations  at 
a  frequency  of  one  cycle  per  day.  (The  spectra  at  the  lowest 
frequency  of  the  low  frequency  segments  of  the  range  and 
refractivity  spectra  may  not  be  correct) .  The  area  vmder 
the  Maui  range  spectrum  is  of  the  order  of  Im.  ,  since  the 
range  changed  less  than  1.3m.  at  Maui  during  a  five  day 
period.  This  range  error  is  less  than  the  range  error  for  a 
spherically  stratified  atmosphere. 

Range  Rate  Errors 

The  range  rate  errors  were  computed  as  a  simple 
/  average  for  a  time  interval  t©  according  to  the  formula 
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(2-17) 


^  ^  R(t  +  O  -  R(t) 

ave 

Table  2-3  contains  RMS  values  of  the  above  quantity  for  periods 
of  time  long  compared  to  to  . 

Azimuthal  Angular  Errors 


The  azimuthal  angular  errors  were  computed  from  the 
expression 


d  01  -  ^  02 
2^ 


X. 

B  cos  P 


where  A  and  A  02  are  the  phase  errors  recorded  at  two 
antennas  separated  by  a  distance  B ;  cos  |3  =  1 ;  and  \  is  the 
wavelength  of  the  transmitted  radio  wave.  RMS  values  of  the 
above  quantity  were  computed  from  instamtaneous  values  of 
the  phase  errors ;  that  is ,  no  time  smoothing  was  done ;  and 
the  RMS  values  were  computed  for  periods  of  data  of  15  min. 
to  1  hour  .  The  RMS  errors  are  plotted  as  a  function  of  the 
baseline  length  in  Figure  2-7.  (Please  excuse  the  introduc¬ 
tion  of  the  English  measure  of  length  on  the  following  figures. 
Sufficient  time  was  not  available  to  put  the  data  in  metric 
units  and  redraw  the  figures  for  this  report.)  The  azimuthal 
errors  decrease  monotonically  with  increasing  baseline  length. 
Note  that  the  results  are  essentially  independent  of  A.  ,  and 
may  be  described  in  terms  of  a  '^range-difference  error" 
divided  by  the  baseline  B.  The  maximum  value  of  the  error 
^Ar  ^PPi^o^ches  0.1  ft  as  B  approaches  10,000  ft. 

Linear  smoothing  of  the  azimuthal  errors  for  the 
longer  baselines  (B  ^  1000  ft.)  is  ineffective  for  a  smooth¬ 
ing  time  tg  less  than  one  minute.  On  the  other  hemd,  the 
phase  difference  power  spectra  curves  for  the  short  baseline 
(B  =  2.2  ft.)  suggests  that  smoothing  for  20  sec.  to  100  sec. 
could  be  very  effective  in  reducing  the  azimuthal  error. 

Azimuthal  Angular  Rate  Errors 

The  azimuthal  angular  rate  error  has  been  calculated 
from  the  azimuthal  error  A^  according  to  the  expression 

4^(t  +  t  )  -  A  pit) 

- rs -  !  (2-18) 
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s^oud  amount 

t  in  aec/\v^ 

O 

Small 

Average 

Large 

5 

.015 

.02 

.055 

20 

.01 

.015 

.03 

60 

.003 

.009 

.02 

- 

TABLE  2-3  The  range  rate  errors  in  centimeters 
per  second  versus  smoothing  time  and 
cloud  zunounts 


Separation  JBaseline,  B  (ft) 

Figure  Z"  7  .  Measured  Variation  of  Azimuthal  Angular  Position  Error 
as  a  Function  of  Interferometer  Baseline  JLength. 


that  is,  the  error  is  an  average  for  the  time  interval  t^. 
The  root  mean  square  of  this  statistic  is  plotted  on  Figure 
2-8  for  a  smoothing  time  t^  =  20  sec.  The  "range-rate 
difference  error"  0/^r  is  also  indicated  by  reference  to  the 
diagonal  dashed  lines.  The  maximum  value  of  0  is  about 
0.001  ft. /sec  for  this  data. 

Table  2-4  is  presented  to  show  that  the  rate  errors  decrease 
as  the  smoothing  time  (t^)  increases . 


TABLE  2-4.  Median  azimuthal  angular  rate  errors  in 
microradians  per  sec  as  a  function  of  baseline  length 
(B)  and  smoothing  time  (t^) . 


N.  B  in  feet 

bQ  in  sec>v^^ 

2.2 

1914 

4847 

5 

19 

0.7 

0.25 

20 

5.5 

0.45 

0.15 

60 

1.7 

0.18 

0.13 

Elevation  Errors 


If  the  refractivity  of  the  atmosphere  is  sufficiently 
isotropic,  then  the  elevation  angular  position  and  rate  errors 
can  be  deduced  from  the  measured  values  of  the  azimuthal  errors 
In  this  case,  the  elevation  errors  are  determined  from  the 
azimuthal  error  data  by  using  the  effective  baseline  length 
instead  of  the  distance  between  the  antennas. 

2.  Cumulus  Clouds 

The  Maui  experiments  did  not  measure  the  propagation 
errors  caused  by  large,  towering  cumuli.  R.  M.  Cunningheun 
of  the  Air  Force  Cambridge  Research  Center  has  made  an  exten¬ 
sive  study  of  the  propagation  errors  caused  by  cumuli.  He 
has  measured  the  refractivity  in  and  arovind  clouds  with  an 
aircraft.  With  these  data  and  cloud  photographs  he  has 
0^  constructed  vertical  cross  sections  of  the  refractivity. 

Then  he  has  computed  the  range  errors  for  various  propagation 
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Separation  Baseline,  b  (ft) 

Figure  2—8  .  Measured  Variation  of  Aximuthal  Angular  Rate  Error 
as  a  Function  of  Interferometer  Baseline  Length, 
20-Sec  Smoothing. 


pathes  along  the  cross  section. 

The  range  errors  which  are  caused  only  by  cumulus 
clouds  are  shown  on  Figure  2-9.  The  cloud  amounts  referred 
to  in  the  caption  show  the  fraction  of  the  earth's  surface 
that  would  be  shadowed  by  clouds  if  the  sun  were  at  the  zenith. 
The  range  errors  are  largest  when  the  clouds  are  largest. 

The  range  error  decreases  as  the  path  length  through  the  clouds 
decreases,  or  as  the  elevation  angle  increases  (theory  predicts 
the  variation  proportional  to  Vcscfi,  plotted  in  Pig.  2-9) . 

The  standard  deviation  of  the  range  error  measured  at  Maui 
during  a  period  of  3^  days  is  also  plotted  and  is  about  the 
Scune  value  as  the  error  caused  by  clouds. 

Cunningheun  computed  rate  errors  from  his  cross  sectional 
data  under  the  assumption  that  the  clouds  and  associated 
refractivity  were  "frozen"  emd  moved  with  the  mean  air  flow, 
which  was  only  in  the  direction  of  the  cross  section.  The 
range  rate  errors  caused  by  cumulus  clouds  are  shown  in  Table 
2-5. 

TABLE  2-5.  The  stamdard  deviation  of  the  range  rate 
errors  in  cm  per  sec  caused  by  cumulus  clouds.  The 
elevation  angle  of  the  propagation  path  is  7°. 


As  before,  t^  refers  to  an  effective  smoothing  time  that  the 
data  were  filtered  linearly.  The  horizontal  speed  of  the 
clouds  in  the  direction  of  the  cross  section  is  given  by  u. 
The  errors  caused  by  cumuli  have  about  the  same  values  as  the 
Maui  errors,  which  were  given  in  Table  2-3. 

Cunningham  computed  the  azimuthal  angular  errors  for 
an  interferometer  from  (2-7),  where  sin  |3'  s  i.  The  standard 
deviation  of  the  errors  are  plotted  as  a  function  of  the 
distance  between  two  interferometer  antennas  in  Figure  2-10. 
The  errors  are  not  smoothed  with  respect  to  time.  The  errors 
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decrease  with  increasing  baseline  length  and  are  seen  to  be 
larger  than  the  maximum  RMS  Maui  errors,  approaching =0.5 
foot  for  the  longest  baseline  with  large  cumuli. 

The  azimuthal  angular  rate  errors  which  are  caused  by 
cumuli  are  shown  for  a  2000  ft.  baseline  interferometer  system 
in  Figure  2-11.  The  errors  show  no  clear  dependence  on  ele¬ 
vation  angle.  The  cumuli  errors  are  larger  than  the  average 
RMS  Maui  errors.  Cunningham's  computations  also  showed  that 
the  angular  rate  errors  caused  by  cumuli  were  inversely  pro¬ 
portional  to  the  baseline  length,  when  that  length  was  between 
2000  ft.  and  25,000  ft. 

Cunningheun  is  presently  computing  the  propagation 
errors  from  refractivity  measurements  made  when  no  clouds 
were  present  over  the  Atlantic  Missile  Rzmge.  When  no  clouds 
were  present,  he  has  measured  horizontal  refractivity  (N) 
changes  as  large  as  60  M-units  in  a  few  meters .  The  changes 
appear  to  be  confined  to  a  layer  of  the  order  of  100  m  thick 
at  a  height  of  about  1  km.  to  2  km.  above  sea  level. 

During  the  previous  discussions  of  fluctuating  propa¬ 
gation  errors,  the  range  rate  and  angular  rate  errors  have 
been  given  for  the  case  that  the  radio  system  is  tracking  or 
guiding  a  stationary  target.  The  rate  errors  also  depend 
on  the  speed  that  the  radio  rays  pass  through  the  atmosphere 
while  following  a  moving  target.  These  errors  are  signifi¬ 
cant  when  the  rays  are  cutting  through  the  atmosphere  at 
speeds  above  that  of  the  prevailing  wind  (which  was  about 
10  ft/sec  for  the  data  shown  above) .  Further  discussion  of 
this  appears  in  Section  4.  of  this  report . 

2.3  Effects  of  the  Ionosphere 

"Radio  seeing"  is  a  problem  of  considerable  interest 
which  affects  many  important  applications .  Just  as  in 
"optical  seeing, "  the  radio  source  near  the  horizon  twinkles 
as  a  result  of  amplitude  scintillation  of  the  arriving  waves, 
and  in  addition,  the  apparent  direction  of  arrival  jitters 
around.  It  is  to  be  expected  that  the  radio  waves  arriving 
from  directions  near  the  horizon  will  be  drastically  modified 
by  the  horizontally  stratified  atmosphere.  There  will,  of 
course,  be  an  overall  gross  bending  or  refraction  of  radio 
waves,  depending  on  the  angle  at  which  they  enter  the  earth's 
atmosphere.  The  purpose  of  this  discussion  is  to  evaluate 
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2^  Florida,  1957j  Large  cumuli;  amotuit  =  0.37^;  t^  =  33b. 

0  Florida,  1957;  Small  cumuli;  amount  =  0.l4;  t^  =  33s. 
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0 

t^  =  408. 

Figure  2>11.  Azimuthal  RMS  angular  rate  errors  as  a  function  of  elevation 
angle.  2000  ft.  hasellne  Interferoneter  system. 
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how  the  distortion  of  radio  fields  by  the  ionosphere  affect 
radio  frequency  positioning  systems,  and  to  specify  the  limits 
on  angle  and  range  corrections  possible  for  ionospheric  re¬ 
fraction  . 

If  the  energy  follows  the  radio  ray,  then  the  refraction 
effects  of  the  troposphere  and  ionosphere  are  additive.  In 
the  troposphere  the  index  of  refraction  is  independent  of 
radio  frequency  and  the  ray  bending  is  a  monotomically 
decreasing  function  of  elevation  angle.  In  the  case  of  iono¬ 
spheric  bending,  the  magnitude  of  the  refraction  is  frequency- 
dependent  and  it  is  a  monotonically  decreasing  function  of 
elevation  angle  only  for  sources  at  distances  large  compared 
to  ionospheric  heights.  For  vehicle  heights  between  one 
hundred  and  four  thousemd  miles,  the  refraction  error  initially 
increases  with  the  elevation  angle,  attains  a  maucimum  value 
at  elevation  angles  on  the  order  of  100  to  200  milliradians, 
and  then  gradually  decreases.  At  ionospheric  heights,  the 
elevation  angle  at  which  the  ionospheric  bending  is  maximum 
is  roughly  proportional  to  the  square  root  of  the  layer 
height.  In  order  to  evaluate  the  remge  and  angle  errors  to 
be  expected,  a  model  of  the  ionosphere  must  be  selected. 

a.  Refractive  Index  of  the  Ionosphere 

As  a  first  approximation,  it  is  assumed  that  the 
index  of  refraction  in  the  ionosphere  is  spherically  strati¬ 
fied  with  a  radial  gradient  only.  Weisbrod  (1959)  has  de¬ 
veloped  a  simple  method  for  computing  atmospheric  refraction 
effects  on  radio  waves,  where  the  refraction  index  profile 
has  been  replaced  by  a  finite  number  of  linear  segments 
whose  thickness  is  small  compared  with  the  earth's  radius. 

Since  there  is  no  limitation  on  the  geometry  of  the  calcu¬ 
lations  related  to  the  shape  of  the  refractive  index  profile, 
the  method  has  a  wide  application  to  other  refractive  effects 
such  as  retardation,  doppler  error  and  Faraday  rotation.  The 
earth's  magnetic  field  has  a  negligible  effect  on  the  phase 
velocity  of  very  high  frequency  radio  waves ;  on  the  other 
hand,  the  Faraday  rotation  may  be  appreciable. 

The  relationship  between  the  index  of  refraction,  the 
radio  frequency  and  the  electron  density  in  the  ionosphere 
is  the  following:  ^ 

n  =  [l  -  ,  (2-19) 
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where : 

N  =  electrons  per  cvibic  meter, 
e  =  electronic  charge  (1.60  x  10"  ^  coulomb), 
m  =  electronic  mass  (9.08  x  10"  kilogreun)  , 

0)  =  2jr  times  the  frequency, 

=  permittivity  of  free  space  (8.854  x  10"^^ 

farad/meter) . 

The  electron  density  variation  with  elevation  is  not  the  same 
all  over  the  earth;  it  varies  with  latitude,  season,  time  of 
day,  sunspot  cycle,  2md  sporatic  conditions.  A  great  deal  of 
information  is  available  on  the  electron  density  profiles 
below  the  height  of  meocimum  density;  these  data  are  computed 
from  ionospheric  vertical  soundings  taken  at  a  network  of 
stations.  Jones  (1962)  has  attempted  to  represent  the  diurnal 
and  geographic  variations  of  these  data  by  numerical  methods. 

He  discussed  the  various  problems  associated  with  the  basic  data 
emd  the  analysis.  To  gain  information  about  the  electron 
profile  above  the  F  maximuAi  requires  a  different  kind  of 
measurement.  Currently,  rocket  and  backscatter  data  are 
available;  in  the  near  future  a  satellite  carrying  a  "top¬ 
side"  sounder  will  yield  ionogrcun  data  above  the  F  meucimum. 

An  electron  density  profile  obtained  from  rocket 
measurement  above  Wallops  Island,  Virginia  at  a  latitude  of 
38^  made  during  the  daytime  up  to  an  altitude  of  600  kilometers 
is  shown  in  Figure  2-13  (taken  from  Jackson  (1961)).  Bowles 
(1962)  has  reported  results  of  observations  of  incoherent 
scatter  obtained  at  the  Lima  Radar  Observatory,  located  at 
minus  twelve  degrees  latitude.  Two  electron  density  profiles 
taken  at  this  station  are  shown  in  Figure  2-14.  These  results 
indicate  a  tendency  of  the  decay  of  electron  density  with 
height  of  the  top-side  of  the  F  region  to  be  exponential.  The 
slope  of  the  exponential  decay  appears  to  vary  from  one  day 
to  the  next  and  from  one  time  of  day  to  the  next  by  a  factor 
of  two.  Electron  density  profiles  available  above  the  F  layer 
maximum  indicate  that  a  hyperbolic  secant  is  a  better  fit  than 
the  Chapman  distribution,  which  is  parabolic.  Therefore,  the 
model  for  electron  density  selected  consists  of  a  parabolic 
variation  below  the  height  of  maximum  electron  density  matched 
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to  a  hyperbolic  secant  profile  above  the  mzuclmvini. 

An  early  study  of  Ionospheric  effects  on  tracking 
(Pflster  and  Keneshea,  1956)  used  simplified  models  of  day 
and  night  profiles ,  Indicated  by  the  rectangular  plots  on 
Fig.  2-14.  The  results  of  this  study,  which  Indicate  the 
order  of  magnitude  of  Ionospheric  range  and  angle  errors, 
are  shown  In  svimmary  form  In  Figs.  2-15  zmd  2-16.  It  should 
be  emphasized,  and  will  be  shown  below,  that  these  errors 
are  subject  to  large  variations  on  a  day-by-day  basis.  In 
addition  to  the  diurnal,  seasonal  and  sunspot-cycle  trends. 

Ionospheric  Model 

For  the  purpose  of  computation  of  refractive  effects 
from  lonosonde  data.  It  Is  necessary  to  postulate  a  model  of 
the  Ionosphere  which  would  approximate  the  observed  data.  The 
normally  available  data  contain  Information  from  which  one  can 
obtain  the  height  of  the  base  of  the  layer,  h^, the  height  of 
the  maximum  electron  density,  h^^,  and  the  critical  frequency 
of  the  layer,  f  .  It  Is  therefore  desirable  to  choose  a  model 
which  has  three^degrees  of  freedom  and  Is  also  In  accord  with 
available  experimental  data  regarding  the  electron  density 
profiles. 


The  shape  of  the  Ionospheric  electron  density  profile 
below  the  maximum  can  be  fairly  well  approximated  by  a  para¬ 
bolic  distribution.  The  shape  of  the  profile  eddove  the  region 
of  the  maximum  density  Is  less  well  known.  From  data  obtained 
from  rocket  soundings  of  electron  density.  It  Is  believed 
that  the  electron  density  above  the  peak  of  the  F  region  does 
not  fall  off  as  rapidly  as  It  might  have  been  expected  from 
the  Chapman  distribution.  Also,  the  Faraday  rotation  experi¬ 
ments  Indicate  that  the  total  electron  content  above  the 
meuclmvun  density  Is  about  three  times  as  large  as  below  It. 
Using  these  facts  the  following  model  may  be  postulated. 

N^/Nq  =  1  -  (1  _  <r)^  ,  0<  <  1  (2-21) 

=  sech  1/47T  (0-1)  ,  >  1 


where : 

Ng  =  electron  density  per  cubic  meter, 
Nq  =  meuclmum  density, 

^  “  (b  —  ho)  /ym  , 
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4; 


yjjj  =  half  thickness  of  the  parabolic  layer 
“  ^  "  ^o  , 

h  s  height  above  the  ground; 
h^  =  height  of  the  base  of  the  layer, 
hj^  =  height  of  the  maximum  electron  density. 

This  model  has  the  following  desirable  characteristics : 

1.  The  model  has  three  degrees  of  freedom,  (h^,  yj^, 
and  Nq  )  which  can  be  obtained  from  ionogramic 
data.  These  pareuneters  vmiquely  specify  the 
entire  distribution. 

2.  The  distribution  is  parabolic  below  the  maxlmvim 
density,  nearly  twice  as  thick  parabolic  immediate¬ 
ly  above  the  maximum  and  exponential  at  great 
heights . 

3 .  The  electron  content  of  the  distribution  above 
the  maoeimum  is  three  times  that  below  it. 

4.  The  entire  electron  density  profile  and  its  de¬ 
rivative  are  continuous  everywhere. 

Figure  3-17  is  a  plot  of  the  ionospheric  model.  The 
heights  of  the  base,  the  maximum  density. and  the  point  of 
interest  define  v ,  and  the  ionospheric  unit  is  obtained 

from  .  2  c 

=  1/2  (N  /N  )  (f  /f)  X  10^  ,  (2-22) 

WllGlTG  * 

f  =  critical  frequency  of  the  layer 
c  1/2  6 

=  8.97  No  '  X  10”  megacycles  per  second, 
f  =  signal  frequency  in  megacycles  per  second. 

The  hg,  hjn»  and  f^,  parameters  refer  to  the  P  layer. 

Using  this  model,  the  refractive  effects  of  the  D 
and  E  layer  are  not  singled  out.  The  reason  for  this  is 
that  they  are  quite  small  in  comparison  with  those  due  to  the 
F  layer  and  are  approximately  accounted  for  by  allowing  the 
electron  density  at  the  bottom  edge  of  the  P  layer  to  be  zero. 
Furthermore,  the  shape  of  the  electron  density  profile  above 
the  mcucimum  is  not  too  well  known  and  since  this  region, 
as  far  as  the  refractive  effects  of  the  ionosphere  are  con¬ 
cerned,  is  probably  much  more  important  than  the  D  and  E 
layers,  it  was  felt  that  the  introduction  of  a  more  complicated 
ionospheric  model  is  not  justifiable.  However,  if  and  when 
accurate  ionospheric  data  from  rocket  soundings  are  available 
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with  comparable  regularity  to  present  radio  soundings,  it  is 
possible  that  the  use  of  more  intricate  models  should  be 
undertaken . 


b.  Computation  of  Ray  Bending  Due  to  the  Ionosphere 


Computational  methods  have  been  developed  which 
afford  simple  means  of  calculating  ionospheric  bending  from 
electron  density  profile  data.  These  are  described  by 
Weisbrod  (1959) .  In  essence  these  methods  are  based  on  the 
following : 

1.  Index  of  refraction  profiles  are  computed  from 
electron  density  profile  data,  where  the  index  of  refraction, 
”n",  as  a  function  of  transmission  frequency,  "f",  and  elec¬ 
tron  density,  "N^,  is  given  by  equation  2-19  and  may  be 
approximated  by 

40. 3N 

n  «  1 - -  ,  (2-23) 

f2 


for  frequencies  well  above  the  critical  frequency.  If  "N" 
units  are  employed  in  the  calculations  where  n  =  1+N  x  10“°, 
then 

-5 

N  =  -  4.03  1-^' 

^2 


X  10 


(2-24) 


In  these  calculations,  the  magnetic  field  of  the 
earth  is  neglected,  since,  at  frequencies  above  100  me,  the 
effect  of  the  terrestrial  field  on  refractive  bending  is 
negligible . 


2.  To  afford  practical  means  of  computation,  index 
of  refraction  profiles  are  approximated  by  straight  line 
segments . 


3 .  Account  is  taken  of  the  spherical  stratification 
of  the  ionosphere.  Thus  the  angles  of  importance  in  appli¬ 
cation  of  Snell's  law  are  those  angles  between  the  ray  tangent 
and  layer  tangent  at  a  given  point.  See  Figure  2-18. 

Principle  formulae  in  the  resulting  calculations  for 


bending  are 


2(h  .  -  h.)  X  10^  N  -  N 

3  ■■  SS  I  >''*'''■<>  "  "  A'"  * 

tan  ^  +  tan  p  500  [tan  P  ,  +  tan^  ] 

J  ^  3  (2_25) 
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where  y  is  the  incremental  bending  in  milliradians  through 
altitud^^'change  from  h .  to  h  (k  referring  to  upper  boundary)  , 
and  the  fi'a  are  the  ray  inclinations  at  the  layer  boundaries. 
The  value  of  /J  at  each  layer  boundary  is  determined  from 
Snell's  law, 

np  cos/3  =  h^a  cos  Sq  ,  (2-26) 


where  a  is  the  earth's  radius,  and  p=  a  +  h,  the 
symbol  denoting  the  height  of  the  layer  boundary, 
bending  through  the  ionosphere  is  then  given  by 


m 
=  £ 


N,  -  N. 


k+1 


total  ]c=0  500  ^  tan  ^  ^Ic 


latter 

Total 


(2-27) 


For  ionospheric  bending,  the  minimum  value  of  /3  even 
for  a  tangentially  departing  ray  (6©  =  0)  is  about  200  milli¬ 
radians.  Under  these  conditions  the  difference  between  the 
ray  inclination  angles  of  the  refracted  and  the  unrefracted 
rays  is  very  small  such  that  very  little  error  is  introduced 
if  equation  (6)  is  written  as. 


m 


\  -  Vi 


=  1,  - 

k=0  500 

tan  0,  , 

+  tan  0, 

f 

• 

k+1 

k 

(2.21?) 


where  the  0's  denote  ray  inclination  angles  at  the  boundaries 
of  the  layer  neglecting  curvature  of  the  ray  within  the  layer. 

Computation  of  Error  Angle 


As  will  be  noted  from  Figure  2-18,  the  apparent  shift 
in  position  of  the  target  is  not  described  by  the  ray  bending 
y,  but  by  the  angle  Moreover,  since  d  is  the  practical 
angle  of  interest  at  this  point  in  the  discussion,  it  is 
necessary  to  consider  the  troposphere  effects  together  with 
those  produced  by  the  ionosphere. 

In  order  to  determine  the  value  of  d  from  calculated 
values  for  y,  it  is  useful  to  express  the  refractive  bending 
in  terms  of  the  angle  subtended  at  the  earth's  center  between 
the  refracted  and  the  unrefracted  rays.  This  angle,  e,  is 
given  by 


39 


(2-29) 


€=  i  -  {e  -  fi) 

=  «  -  (Nq  -  N)  cot  0 


where : 


€  = 

6  = 

= 

N  = 


- 1  £  i 

surface  value  of  the  refractivity, 

value  of  the  refractivity  at  the  target  height. 


At  infinite  distances  cot  6  approaches  zero  and  c 
and  6  become  equal  to  each  other.  For  computational  reasons 
it  is  usually  convenient  to  split  equation  (2-27)  into  the 
tropospheric  and  ionospheric  components. 

-  Nq  cot  e  (2-30) 


=  5^  +  cot  e 


(2-31) 


The  quantity  of  the  greatest  practical  interest  is  the 
elevation  angle  error  6  ,  which  can  conveniently  be  expressea 
in  terms  of 


€  tan  e  t  €  ^/2 

i~+'TanTrTrTwr“0r 

o 


(2-32) 


It  should  be  noted  that  the  tropospheric  and  the 
ionospheric  contributions  to  5  are  not  strictly  additive. 
However,  in  nearly  all  practical  cases  c/2  is  much  less 
than  tan  6  and  c  is  much  less  than  tan  6  -  tan  6o,  so  that 
only  a  negligible  error  is  introduced  if  0  is  considered  to 
be  directly  proportional  to  c .  This  approximation  is  ex¬ 
tremely  convenient  since  it  permits  a  separate  treatment  of 
the  tropospheric  and  the  ionospheric  5's. 


Range  Errors  and  Faraday  Rotation 

The  range  errors  can  be  evaluated  by  a  method  analo¬ 
gous  to  that  used  for  computing  y: 


Ar^  =  t 


m 


'-K  *  Vi> 


-  V 


v_n  1000  (sin  B  +  sin  0  ) 
k+1  k 


(2-33) 


40 


where : 


h  =  height  in  kilometers, 
and  4r^  =  range  error  in  meters. 

The  accuracy  of  equation  (2-33)  is  adequate  for 
most  practical  purposes .  For  layer  leuninations  of  less 
than  100  and  a  tangentially  departing  ray,  the  values 
of  Ar^  are  within  5%  of  the  exact  values.  For  thinner 
Izuninations ,  (which  is  normally  the  case)  or  higher  angles 
of  elevation,  the  agreement  is  even  better. 

Equation  (2-33)  also  offers  a  convenient  method  for 
computing  the  Faraday  rotation.  It  turns  out  that  the 
number  of  rotations  of  the _ plane  of  polarization  is  very 
nearly  proportional  to  Ar^ .  For  the  case  of  a  thin  ].ayer 
the  relationship  between  the  Faraday  rotation  and  A is 


aq 


2  cos  <f> 


(2-34) 


where : 

O  =  the  number  of  rotations  of  the  plane  of  polari¬ 
zation  for  a  double  passage  through  the  layer, 

4^  =  the  angle  between  the  wave  normal  and  the  magnetic 
field, 

=  gyro  wavelength. 

Both  0  and  Xg  are  functions  of  a  position.  However, 
the  rate  of  change  of  these  quantities  is  sufficiently  slow 
so  that  F  region  values  of  these  parameters  may  be  treated 
as  constants  applicable  to  the  entire  path. 

c.  Statistics  of  Ionospheric  Errors 

To  illustrate  the  magnitude  and  variation  of  the 
errors  in  elevation  angle  and  range  which  may  be  ascribed 
to  the  ionosphere,  the  results  of  calculations  are  given  in 
Figures  2-19  through  2-22  for  100  me  for  three  different 
locations:  Iverness  or  Leuchars,  Scotland;  Thule,  Greenland; 
and  Fairbanks,  Alaska.  Since  ionospheric  parameters  are 
subject  to  wide  fluctuations,  it  is  necessary  to  select 
representative  situations  and  group  them  in  some  logical 
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Flgun  2-21.  Diatribution  of  ionoapharic  range  error  Ar”^ 

(fsioo  mcpa,  6q*0,  h»1000  n.od.) 

Low  aunapot  activity. 


fashion.  The  three  moat  important  parameters  are  the  sunspot 
cycle,  the  seasonal,  and  the  diurnal  variations.  Thus,  data 
is  shown  for  March,  June,  September,  and  December  of  1954  and 
1957.  Individual  points  are  for  particular  days  of  the  month. 
Average  diurnal  variations  may  therefore  be  inferred  from 
maximum  point  density  distribution.  The  years  of  1954  and 
1957  represent,  respectively,  low  and  high  sunspot  activity. 
Since  there  is  no  1954  ionospheric  data  for  Thule,  the  data 
for  Resolute  Bay,  Canada  were  used  in  the  calculations  for 
that  year  since  the  location  of  Thule  and  Resolute  Bay  lie 
in  reasonable  comparable  geographic  and  geomagnetic  latitudes 
and  thus  ionospheric  data  from  the  two  locations  are  similar. 
Since  the  data  in  Figures  2-19  through  2-22  were  calculated 
for  100  me ,  to  employ  these  results  for  higher  transmission 
frequencies,  advantage  may  be  taken  of  the  inverse  square 
dependence  of  ionospheric  index  of  refraction  on  frequency. 
Thus ,  error  angles  for  200  me  could  be  determined  from  these 
data  by  scaling  down  each  value  for  error  angle  by  a  factor 
of  four . 


Figures  2-23  and  2-24  show  the  change  in  error  angle 
and  range  error  with  sighting  elevation  angle,  Oq.  These  data 
were  computed  for  a  critical  frequency  of  one  megacycle.  Thus, 
the  values  f or  and^r^in  Figures  2-19  through  2-22  which  are 

given  for  9q  =  0,  may  be  determined  for  other  sighting  eleva¬ 
tion  angles .  It  is  interesting  to  note  the  increase  of  ^ 
with  Oq  to  a  meucimum  near  100-200  milliradians . 

The  following  interesting  generalizations  may  be 
drawn  from  the  data  in  Figures  2-19  through  2-22: 

1.  When  there  is  a  diurnal  variation  it  is  strongest 
in  December  and  weakest  in  June. 

2.  The  shape  of  the  diurnal  variation  is  the  same 
throughout  the  sunspot  cycle,  but  the  values  of 
the  errors  are  approximately  five  times  as  great 
during  the  sunspot  maximum  as  during  the  sunspot 
minimum. 

3.  The  diurnal  variation  is  greatest  at  the  station 
farthest  from  the  geomagnetic  pole  and  almost 
non-existent  near  the  pole. 

4.  The  relative  spread  of  values  at  any  one  hour  is 
greater  during  the  night  time  hours  than  during 
the  daylight  hours  and  is  larger  near  the  pole 
than  away  from  the  pole. 

5.  The  diurnal  variation  may  exceed  an  order  of 
magnitude . 
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Of  the  three  factors  which  affect  the  magnitude  of  and 
^r^fthe  critical  frequency  of  the  layer  appears  to  be  the 
most  dominant,  It  is  for  this  reason  that  the  diurnal  shapes 
of  i  ^  and  ^r^  are  quite  similar.  The  second  factor  which 
influences  the  magnitude  of  the  error  is  h^.  The  lower  the 
layer  the  more  oblique  is  the  angle  of  entrance  into  the 
layer  and  consequently  the  greater  the  error,  nie  third 
factor,  has  a  greater  effect  on  dr^  than 

In  general  the  critical  frequency  has  a  more  pro¬ 
nounced  diurnal  variation  in  winter  thzui  in  summer,  is  greater 
at  Inverness  than  at  Thule,  and  is  generally  higher  during 
the  peak  of  the  sunspot  cycle  tl^an  at  its  minimum.  However, 
the  effects  of  h^  on  6^  and  d r^  cannot  be  neglected.  In 
December,  h^  decreases  during  the  day  and  Increases  at  night 
approximately  in  phase  with  f^.  The  condoination  of  effects 
greatly  accentuates  the  diurnal  variation  in  and  dr^. 

In  June,  h^  is  approximately  constant,  but  the  variation 
which  does  occur  is  in  such  a  direction  as  to  partially 
cancel  out  the  smaller  effects  of  f^,  resulting  in  almost 
no  diurnal  variation  of  The  variations  in  y,Q  have  a 

similar  effect  on  ^r^. 

In  addition  to  the  errors  discussed  above,  short¬ 
term  fluctuations  in  electron  density  will  cause  variation 
in  range  measurements  over  periods  of  ten  seconds  to  a  few 
minutes.  Figure  2-25  is  em  estimate  of  such  fluctuations, 
based  on  data  from  Pfister  and  Keneshea  (1956)  euid  Colin 
(1960)  . 


The  smooth  gross  structure  of  the  N  profile  causes 
errors  which  can  be  corrected  if  one  knows  the  shape  of  the 
N  profile,  but  the  inhomogeneities  cause  an  Indeterminancy 
in  the  actual  vehicle  location  which  cannot  be  resolved. 

A  discussion  of  observed  elevation  angle  errors  using  a  sea 
inter fercnne ter  technique  is  given  in  Appendix  A. 


f 
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fluctuation  In  feet 


Figure  Z-25.  Ionospheric  range  fluctuation  vs.  Frequency 
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3.  ATMOSPHERIC  SOUNDING  FOR  CORRECTION  OF  TRACKING  DATA 
3.1  Methods  of  Measuring  &  Estimating  Tropospheric  Profiles 

The  index  of  refraction  profile  may  be  obtained  indirectly 
through  use  of  radiosonde  or  dropsonde  instruments .  Although 
there  are  a  nundser  of  classes  of  these  instruments,  their 
function  is  similar.  Carried  aloft  on  free  balloons,  or 
jettisoned  by  parachute  from  aircraft  or  rochets,  the  device 
saunples  the  absolute  values  of  dry-bulb  temperature,  relative 
humidity,  and  in  some  instances,  atmospheric  pressure  of  the 
air  through  which  they  move  in  generally  oblique  directions 
as  a  result  of  lateral  displacement  by  winds .  The  measure¬ 
ments  are  transmitted  as  a  modulation  of  an  r-f  carrier 
signal.  A  switching  device  alternately  switches  circuits 
to  the  transmitter,  thereby  providing  a  sequence  of  samples 
of  the  parameters  being  measured,  one  at  a  time.  It  is  thus 
necessary  to  interpolate  over  the  interruptions  in  these 
data  by  arbitrarily  assuming  that  the  intervening  variations 
are  smooth.  The  data  are  normally  recorded  on  a  strip-chart 
(ground  based  or  in  an  aircraft)  as  a  deflecting  trace.  De¬ 
flections  must  be  measured  and  properly  interpreted  to  recon¬ 
struct  the  original  parameters.  There  may  be  an  associated 
radar  tracking  the  sonde  to  give  space  position.  At  many 
installations  electronic  computers  are  utilized  in  determin¬ 
ing  from  the  chart  readings  the  electrical  and  optical  index 
of  refraction,  as  well  as  many  other  atmospheric  parameters. 

In  the  newer  radiosonde  packages  (AN/GMD2) ,  the  tempera¬ 
ture  element  is  a  resistance  type  thermister  (ML-419/AMT-4) . 
The  humidity  element  (ML-476/AMT)  is  carbon  coated.  This 
humidity  element  is  a  vast  improvement  over  the  old  type 
ML-418/AMT  which  was  lithium  chloride  coated.  The  carbon 
element  has  virtually  eliminated  the  lag  in  the  response  to 
humidity  changes  that  was  characteristic  of  the  lithium 
chloride  element . 

Because  of  the  alternate  observations  of  the  parameters , 
and  some  lag  in  the  sensors,  and  since  ducting  layers  may  be 
only  a  few  hundred  feet  in  thickness ,  the  radiosonde  may 
fail  to  detect  the  existence  of  such  conditions,  or  may  simply 
integrate  over  the  whole  range  of  values  between  the  bottom 
and  top  of  a  duct,  so  that  an  indication  of  the  true  gradients 
throughout  the  layer  cannot  be  achieved. 
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Once  a  determination  is  made  of  the  necessary  quantities, 
the  index  of  refraction,  n,  or  more  conveniently  refractivity, 
N,  can  be  determined  from  the  equation 

N»  (n-1)  10^  »  2^  (P4  i|i2E)  (3-1) 

where  T  is  absolute  temperature  in  degrees  Kelvin,  P  is 
atmospheric  pressure  in  millibars,  and  p  is  vapor  pressure. 

The  constants  are  reasonably  correct  as  determined  by  various 
empirical  studies.  These  particular  constants  are  a  re-expres¬ 
sion  of  those  recommended  by  Smith  and  Weintraub  to  yield  an 
over-all  standard  error  of  in  N,  assuming  P,  T,  and  p  to 

be  error  free. 

Much  more  reliable  refraction  data  can  be  obtained 
through  the  use  of  any  one  of  several  types  of  apparatus 
known  as  microwave  refractcxneters ,  which  measure  directly  the 
refractivity  of  the  atmosphere.  Since  the  velocity  of  propa¬ 
gation  of  radio  waves  is  a  direct  function  of  the  radio  re¬ 
fractive  index,  the  resonant  frequency  of  cavity  resonators 
and  tuned  resonant  frequency  devices  in  general,  are  also  a 
direct  function  of  the  radio  refractive  index  of  the  material 
within  the  resonant  cavity  or  the  dielectric  material  in  the 
capacitor  determining  the  frequency  of  resonance.  Refracto- 
meters  using  this  principle  can  be  classed  into  two  broad 
categories :  those  that  are  capable  of  determining  only  the 
relative  change  in  refractivity  from  one  point  to  another  in 
space,  and  those  that  determine  the  absolute  value  of  N  at 
any  such  point.  The  refractometer  can  ascertain  within  plus- 
or-mlnus  a  few  N-unlts  the  small  scale  variations  of  refrac¬ 
tivity  over  the  path  that  it  follows.  The  lag  of  a  refracto¬ 
meter  is  dictated  by  the  speed  of  its  carrier,  and  the  time 
required  to  flush  the  sample  air  from  the  instrument.  These 
usually  combine  to  permit  accurate  sensing  of  layers  as  thin 
as  100  feet. 

AFSC  is  currently  attempting  to  develop  an  operational 
expendable  refractometer,  which  will  make  economically  feas¬ 
ible  the  routine  acquisition  of  highly  accurate  refraction 
data. 
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In  instances  when  it  is  not  possible  to  sample  the 
atmosphere,  approximations  of  various  sorts  may  be  used  to 
estimate  what  the  N-structure  is  but  with  considerable  loss 
in  the  probable  reliability  of  the  profile  and  of  any  cor¬ 
rections  of  space  positioning  data  based  on  the  generated 
N-structure .  Some  of  the  approximations  include  use  of  a 
standard  atmosphere,  a  local  average  atmosphere,  or  one  of 
the  above  adjusted  for  seasonal  and  diurnal  average  varia¬ 
tion.  Use  of  a  suitable  mathematical  model  which  is  a 
function  of  surface  refractivity  (which  Is  readily  determined 
from  surface  observations)  is  perhaps  the  best  approximation 
procedure  to  use. 


i 

1 
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3.2  Ray  Tracing  Prodecures  for  Correction  of  Tracking  Data 

Once  a  table  of  index  of  refraction  versus  altitude  is 
obtained,  one  might  think  it  possible  to  obtain  an  exact 
propagation  correction  for  an  observed  elevation  angle  and 
range.  However,  there  are  several  difficulties.  Refractiv- 
ity  profiles  cannot  be  obtained  for  all  places  and  times  of 
interest.  Assumptions  are  therefore  normally  made  that  the 
atmosphere  is  non-moving,  spherically  stratified,  and  the 
strata  are  homogeneous,  any  or  all  of  which  may  not  be  true. 
Thus,  any  error  in  the  index  of  refraction  profile  resulting 
from  measurement  or  from  the  above  assumptions  will  be  prop>a- 
gated  through  the  calculations,  resulting  in  errors  in  the 
computed  corrections . 

The  equations  used  to  perform  the  calculations  leading 
to  corrections  for  propagation  effects  are  normally  based 
on  "ray -tracing"  theory  which  assumes  a  ray  of  energy, 
rather  than  a  "wave  front"  travels  between  the  target  and 
the  receiver.  Restrictions  of  ray-tracing  include: 

(1)  The  refractive  index  should  not  change 
appreciably  in  a  wave  length. 

(2)  The  fractional  change  in  the  spacing  between 
neighboring  rays  (initially  parallel)  must  be  small  in  a 
wave  length. 

Conditions  (1)  and  (2)  will  be  violated  if  the  gradient 
of  refractive  index,  dn/dh,  is  very  large.  "Trapping"  can 
occur  whenever  a  layer  of  refractive  index  exists  with  a 
vertical  decrease  of  N  greater  than  157  N-units  per  kilo¬ 
meter.  A  layer  of  this  type  is  a  "duct"  and  the  mode  of 
propagation  through  such  a  layer  is  similar  to  that  of  a 
wave  guide  (Booker  and  Walkinshaw,  1946) .  Condition  (1) 
should  be  satisfied  if 

(dn/dh)  per  km  #3  2\ 

N  <  .002  f^^ 

where  refractivity,  N,  is  defined  as 
N»(n-1)  10® 
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and  f  is  the  carrier  frequency  in  kilocycles  (Bean  and 
ThayeP?  1959) . 

In  addition  to  the  ray  theory  assumptions,  some  error 
is  introduced  into  the  developed  equations  through  the 
spherical  earth  assumption,  and  even  more  by  a  flat  earth 
assumption  when  such  is  used. 

The  correction  equations  necessarily  include  definite 
integrals  which  cannot  be  integrated  explicitly.  It  is 
therefore  necessary  either  to  make  further  simplifying  — 
and  error  producing  —  assumptions  which  modify  the  equations 
in  such  a  way  that  they  can  be  explicitly  integrated,  or  to 
choose  a  numerical  integration  technique  to  perform  the  inte¬ 
grations.  All  niunerical  integration  techniques  are  approxi¬ 
mations,  and  there  are  many  techniques,  each  with  possible 
variations  to  choose  from.  Which  technique  is  "best"  for  a 
given  integral  is  not  easy  to  determine.  It  is  thus  apparent 
that  errors  will  be  introduced  because  of  the  nature  of  the 
required  computations. 

Once  a  basic  technique  is  chosen,  there  remain  still 
other  possibilities  of  error.  Many  mathematical  formulas 
are  such  that  for  certain  critical  spreads  of  values,  to 
prevent  loss  of  accuracy  extended  precision  (retaining  a 
larger  nundser  of  significant  figures  in  the  computation 
than  is  available  in  a  ccmiputer  "word")  is  required  in  the 
computations,  and  at  times  the  formulas  may  completely 
fail.  It  is  thus  sometimes  desirable  if  possible  to  re¬ 
express  the  equations  being  used  —  either  in  critical  areas 
or  over  the  entire  range  —  in  a  mathematically  equivalent 
form  less  susceptible  to  computational  error. 

Even  with  these  limitations  of  ray-tracing  theory  and 
numerical  integration  techniques,  experience  has  proven  they 
are  still  capable  of  providing  much  more  accurate  refraction 
corrections  than  are  methods  based  upon  approximate  curves 
or  assumed  standard  conditions  which  do  not  accurately  reflect 
the  actual  weather  conditions  or  properly  utilize  the  geo¬ 
metrical  considerations. 

Understanding  that  any  technique  is  not  error  free,  it 
is  desirable  to  select  one  which  will  produce  a  minimum 
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error.  Several  good  techniques  will  be  described,  essentially 
those  used  by  NBS,  PMR,  AMR,  and  AP6C.  No  comparison  of 
results  of  the  various  methods  against  an  adequate  standard 
has  been  made,  hence,  there  may  be  a  definite  order  of 
relative  merit. 

The  NBS  development  of  correction  equations  is  given 
first. 

Snell's  law  for  a  spherical  earth,  forming  the  basis 
of  ray  tracing  theory,  may  be  expressed  in  the  forms 

n^  (r+ho)  Cos  6©  *  ni(r+hi)  Cos  Si  ■  Const  (3-3a) 

niri  Cos  Si  =  ni+i  Cos  0^41  *  Const  (3-3b) 

where  r^  =  r+h^^ 

r  =  local  radius  of  the  spherical  earth 
hj^  =  height  above  the  spherical  earth 
In  (3 -3b)  let 

ni+i  -ni  +  An,  »  ri  +  4r,  ei4i  *  Bi  +  40 

where  4n,  4r,  40  are  small  increments.  Choosing  these 
increments  sufficiently  small,  they  may  be  considered  as 
infinitesimals.  Then  s\4)stituting  in  (3-3b)  we  have 

niriCos  0^  «  ni4l  r^+i  Cos  ■  (3-4) 

(nj^+dn)  (r^+dr)  Cos  (B^+dB) 

Expanding,  discarding  products  of  differentials,  and  setting 
sin  d0  *  d0,  cos  d0  «  1 

gives 

n^ri  Cos  ©i®  n^ri  Cos  ©i-nir^  sin  0^  dB+n^  Cos  B^dr+riCos  B^^dn 

(3-5) 


SB 


t 


Subtracting  Cos  0^  from  both  sides,  dividing  by  Cos  0i, 

and  solving  the  resultant  equation  for  tan  0  d0  gives 


tan  0^  d0  * 


n^  dr+r^^  dn 


«i^i 


or 


dn 


d0  »  Tf  Cot  0  + 


dr 


dr  dn 


Cot  0 


(3-6) 


(3-7) 


From  the  infinitesimal  portion  of  the  ray  path  pictured  in 
Fig.  (3-2) ,  it  can  be  seen  that 

rjd^ 


Cot  «i  -  ^  *  X  r^  d0 


Cot  0, 


dr 


Thus , 


T-  Cot  e  .  ^  #  *  ■*♦ 


(3-8) 


and  (3-7)  becomes 


d0  B  %  Cot  0  +  d^ 
n 


(3-9) 


From  Fig.  (3-1) , 

(0  4  90)  +  ♦  +  (90  -  0)  4  (180-  T)  -  360 

o 


or 


0©  4  ^  -  T  »  0 


(3-10) 


Differentiating, 


d0  »  d^  -  dT 


Cranparing  this  with  (3-9) ,  it  is  seen  that 

dr  «  Cot  0 

n 


*  ( 


or 


-  ^ 


Cot  0 


(3-11) 


59 


FOA/£TjR  y  OF  THF 


9 


Figure  3-1 


Figure  3-2 


In  the  numerical  evaluation  of  integral  for  T,  the  height 
increments  may  be  chosen  such  that  the  change  in  n  is  never 
large,  so  that  a  linear  gradient  of  n  may  be  assumed  between 
sample  points,  with  negligible  error.  For  the  case  of  a 
linear  n  gradient  between  hj^  and  h,,  the  integral  can  be 
solved  in  essentially  closed  form  by  assuming  that  over  a 
small  increment  ^  and  Cot  6,  can  be  considered  to  be  con¬ 
stants  expressible  as  s  and  Cot  and  this  can  be  factored 
out  of  the  integral.  5ncorporating  this  concept,  equation 
(3-11)  can  be  restated  as 


3  Cot  5. 
«r=  -  X 

i=l  ”i 


n 


where 


0.  = 

1 


®i  ^  »i 


rn. 

L 

dn  =  - 

j 

L 

Cot  S . 

(n 

n  1 

i-1 

i=l 

1 

1 

n. 

n _ L. 

+  n. 

-1 

i 

n  —  " 

2 

1-1) 

(3-12) 


and  6^  is  obtained  through  a  re-expression  of  Snell's  Law 
in  a  half-angle  sine  form  to  permit  retention  of  significance. 


9  -  2  arc  sin  2 


„  6  h-h  n  -n 

F5Tr  -  -|— 

s 


Cos  6 


o, 

(3-13) 


.]) 


Once  T  and  6  for  the  observation  of  interest  are  estab¬ 
lished,  the  elevation  angle  error,  6  ,  may  be  obtained  from 

n  -n 

— - —  +1  -  Cost  -  sin  T  tan 

6  =  T  -  arc  tan  -  (3-14) 

sinT  -  Cost  tan  6  +  ”s  tan  0 

o  T — 


The  true  elevation  angle  is  then 

0^  =  +  «  (3-15) 

the  angle  at  the  center  of  the  earth  is 

^  =  r  +  6  -  0o  .  (3-16) 
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and  the  true  range  (from  geometry  of  Pig.  (3-1))  is 


(  ^ 


R 


t 


(r-fh)  sin  <i> 
Cos  6^ 


(3-17) 


Since  the  NBS  work  is  all  theoretical  the  height  to 
which  integration  is  performed  is  provided  the  computer  pro- 
graun  as  a  constant.  For  the  correction  of  tracking  radar 
data,  the  height  to  which  integration  is  to  be  performed  can 
be  determined  as  the  height  above  a  spherical  earth  indicated 
by  the  observed  elevation  angle  and  range. 


h=r[(l.(|)2.2  (b)  ,i„  .  l] 


(3-18) 


To  increase  accuracy  of  the  correction  the  solution  can  be 
performed  as  indicated  above,  then  a  better  approximation 
of  height  of  the  target  can  be  obtained  from  (3-18)  using  the 
value  obtained  for  0^  in  place  of  0^.  Using  this  height,  the 

last  increment  of  the  nximerical  integration  and  subsequent 
computation  should  be  repeated. 


The  mathematical  equation  formulation  used  by  APGC  is 
patterned  after  that  used  by  NBS,  except  for  the  evaluation 
of  the  integral  for  T  (Eq.  (3-11) )  which  is  done  through 
usage  of  a  ten  point  Gauss-  Legendre  method  of  mechcunical 
quadrature.  This  method  yields  exact  results  for  the  func¬ 
tion  if  it  can  be  expressed  as  a  polynomial  of  degree  nine¬ 
teen  or  less . 


The  AMR  mathematical  development  is  slightly  different. 
Using  the  equation  (3-3a)  expression  for  Snell's  law,  divid¬ 
ing  both  sides  by  r,  and  solving  for  Cos  9  gives 


Cos  0  = 


Cos  0o 

n  (1+  ll) 
r 


k 

n(l4  ll) 
r 


(3-19) 


letting  k  denote  the  constant  numerator  of  the  fraction. 


f 
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Substituting  this  expression  for  Cos  0  into  (3-8) ,  one  gets 


d0  = 


Cos  6  dh 


k  dh 


(r+h)  V  l-Cos2^ 


(r+h)n(l4^)  (14^)2  _  ^2 


n2  (1+$)^ 


k  dh 


k  dh 


(r4h)  Vn2  (14^)2  _  r(l+^)  Vn^  (14$)“^  -It  (3-20) 


where  dh  =  dr  of  equation  (3-8)  .  Integrating, 

^t 

dh 


</>=“/  (1421)  Vn2  (14^)2  -k2 

^  h_  r  ^ 


(3-21) 


AMR  evaluates  0  using  a  five  point  Gauss -Legendre  numeri¬ 
cal  integration  technique,  which  should  yield  exact  results 
if  the  function  under  the  integral  can  properly  be  expressed 
as  a  polynomial  of  degree  nine  or  less . 

From  Pig.  (3-1),  the  law  of  sines  applied  to  triangle 
ABD  gives 

sin  ^  siJf^  (90  4  6^)  Cos  0^ 


r4h 


r4h 


(3-22) 


From  the  definition  of  Cot  the  length  of  side  AC  is 

and  thus  the  length  of  side  ^  is 


CD  =  (r4h)  - 


Cos 
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Now  in  triangle  BCD,  the  law  of  sines  gives 


sin 

-  sir?- 


sin  (90  4  4)  -  Cos  ^ 


or 


(r+h)  Cos  ^  -  r 
sin  0^ 


(3-23) 


Substituting  this  into  (3-22)  and  solving  for  Cot  6^ 


Cot  0^.  = 


(r+h)  sin  ^ 
(r+h)  Cos  ^  -r 


sin  ^ _ 

Cos  <f>  -  1 

h 

(1+p 


(3-24) 


The  true  elevation  angle  is  computed  using  equation  (3-24) . 

The  same  comment  applies  here  as  stated  immediately  after 
Eg.  (3-18) ,  that  one  "iteration"  will  improve  accuracy.  The 
true  range  is  calculated  as  indicated  in  Bq.  (3-17) . 

The  PMR  development  follows  that  of  AMR  through  Eq. (3-21) . 
At  this  juncture  a  substitution  is  made : 


sin  0 


Cos^  0 


n  (14^) 


and  (3-21)  becomes 


0  = 


Cos  0 
sin  0 


n  dh 


(3-25) 


(3-26) 


The  index  of  refraction  profilie  if  plotted  can  be  broken 
into  a  number  of  segments  of  essentially  constant  gradient. 
This  can  be  expressed  as 


n.  =  n.  ,  +  a. 
1  1-1  1 


(hi  -  h.  _i) 


(3-27) 
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where 


"i  -  "i-1 

“  TT—^TT — 


‘i-l 


is  the  slope  of  the  segment.  Differentiating  (3-27) 


d  h. 


a. 

1 


Now,  differentiating  Eq.  (3-19)  with  respect  to  h,  we  have 


-  sin  0 


de 

dh 


Js 

n 


.  -  (1+ii) 

r  r 


hx-2 


n^(l4~) 


dn 

3H 


or 


H  e,  de  =  [a  (i+J)  4  |1  dh 


Cos  0 


where 


,  dn 


which  can  be  rearranged  in  the  form 


„  n _ r  k  sin  0 

-  lalr+hUn) 


(3-28) 


Substituting  this  into  equation  (3-26) , 


0 


-J 


(r4h) 4n 


de 


(3-29) 


Over  appropriate  intervals,  where  ^  is  approximately 
constant,  the  coefficient 


D  =  n 

a(r4h) 4n 

is  approximately  constant.  There  will  be  a  small,  approxi¬ 
mately  linear  variation  of  this  coefficient  over  the  altitude 
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interval  chosen.  The  mean  value  of  the  coefficient  is  very 
nearly  the  average  of  the  coefficients  at  the  beginning  and 
ending  of  this  altitude  interval .  Thus , 

h  (n^  +  n^) 


D. 

1 


"i 


a. (r+h. ) 4n. 
1  11 


n.  -n 
h 


i-^'i-l  f 

V.K~ 

1  1-1  L 


(nj  4  nj.i) (hj  -  hj.i) 


(n 


i  -"i-i*  b  *  <‘'1  * 


+  (n^+n._^) (h^-h^_^) 


(3-30) 


and  equation  (3-29)  becomes 


(3-31) 


and  0.  is  evaluated  from  equation  (3-19) .  The  height  to 
which^integration  is  performed  may  be  determined  through 
usage  of  equation  (3-18) ,  and  the  true  elevation  angle  from 
equation  (3-24) ,  and  the  true  range  from  equation  (3-17) . 

From  a  computational  accuracy  and  convenience  standpoint, 
it  would  seem  that  an  acceptable  standard  method  based  on 
portions  of  all  of  the  above  methods  could  utilize  equations 
(3-18),  (3-13),  (3-30),  (3-31),  (3-24)  and  (3-17).  That  is, 
the  steps  would  be : 

a.  Compute  the  height  to  which  integration  should  be 
performed. 

h=  r  r  ^14(R)^  4  2  (^  sin  0^  -1 
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b.  For  the  upper  end  of  each  increment  chosen  such 
that  ^  is  essentially  constant  over  the  altitude  interval, 
compute  the  local  elevation  angle. 


/  r4h-  r  .  o  0o  "s-"  1  \  ^ 

0=  2  arc  sin  1  2  sin^  ^  4  ^  Cos  0^  J  J 


c.  Compute  the  “integration  constant"  for  this  interval 


D,= 


("i  +  "i-l'  O'!  ->>1.1) 


^  f  +  (n^+n^-1)  (h^ 


d.  Compute  the  angle  at  the  center  of  the  spherical 
earth  determined  by  the  radar  site  and  the  item  being  tracked 
retaining  also  the  values  of  the  two  components  defined  as 
A  and  B. 


0=  D.  (e._  4  Dj(e.-  e._i)  s  a  4B 

i=l 


e.  Compute  the  first  estimate  of  true  elevation  angle. 


0.  =  Arc  Cot 
t  Cos  <p  - 


d"!) 


f.  Recompute  height  from  (a)  using  0^  in  place  of  0q. 

g.  Recompute  Sj,  Bj  from  b  and  c,  and  compute  a  new  0 
using  the  original  component  A  and  a  newly  computed  component 

B. 

h.  Compute  the  final  estimate  of  true  elevation  using  e. 

i .  Compute  the  true  range .  R^.  =  - 

^  Cos  0^ 
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Correction  of  the  elevation  angle  for  optical  tracking 
devices  can  also  be  done  in  the  above  manner,  utilizing 
the  optical  index  of  refraction. 
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3.3 


Estimated  Accuracy  of  Corrected  Measurements. 


The  accuracy  of  atmospheric  refractivity  measurements 
has  been  discussed  by  McGavin  (1962)  who  states  : 

"Meteorological  sensing  is  limited  mainly  by  the  in¬ 
accuracy  in  measuring  humidity  which  under  ideal  conditions 
appears  to  limit  the  accuracy  to  1  1.0  N.  Gradient  measure¬ 
ments  utilizing  radiosondes  reflect  an  accuracy  no  better 
than  1  3  N  units .  Radio  frequency  refractometers  are  capable 
of  accuracies  as  much  as  an  order  of  magnitude  better  than 
that  achieved  by  meteorological  sensors .  Lightweight  re¬ 
fractometers  have  been  devised  for  balloon-borne  and  dropsonde 
measurements  reflecting  accuracies  inferior  to  the  conventional 
refractometer  but  superior  to  the  radiosonde." 

Since  the  surface  refractivity  is  in  the  order  of  330  N 
units ,  the  range  corrections  resulting  from  use  of  complete 
profiles  with  the  above  accuracies  will  be  in  error  by 
0.396  to  1%.  Angle  corrections,  which  depend  upon  gradients 
of  refractive  index  and  are  more  sensitive  to  variations  in 
the  profile,  will  have  a  larger  percentage  error,  usually 
between  2%  and  4%  of  the  initial  value  of  error.  The 
accuracy  of  the  angle  correction  will  be  better  for  long- 
range  targets  than  for  those  which  lie  within  the  troposphere, 
since  the  elevation  error  approaches  the  total  bending  T  in 
the  long-range  case . 

Bean  and  CahOon  (1957)  showed  by  ray-tracing  calculations 
that  the  residual  error  in  r  could  be  reduced  to  within  50 
uradians  RMS  at  elevation  angles  above  3°  using  surface  re¬ 
fractive  index  only.  This  corresponds  to  an  error  of  196  of 
the  initial  value  at  3^.  Radio  measurements  described  by 
Iliff  and  Holt  (1962)  demonstrated  that  good  accuracy  could 
be  achieved  in  prediction  of  total  bending  from  surface 
refractive  index.  Their  results  showed  residual  errors  of 
150  uradians  at  2°  elevation,  or  2.596  of  the  initial  error 
of  6000  u  rad.  Experience  at  the  test  ranges,  reported  to 
the  Panel,  confirms  the  feasibility  of  correcting  tracker 
data  to  about  50  uradians.  Both  AMR  and  APGC  reported 
residual  bias  errors  in  the  order  of  100  uradians  for 
An/PPS-16  radar  data,  and  a  large  portion  of  this  figure 
must  be  attributed  to  the  radar  itself. 
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The  present  state  of  the  art  in  correction  of  range 
and  elevation  data,  using  combined  radiosonde  and  refracto- 
meter  data  to  derive  accurate  surface  refractivity  and 
profiles,  is  estimated  to  provide  the  accuracy  of  correction 
shovm  in  Table  3-1. 


TABLE  3-1. 


Optim\am  Accuracy  of  Range  and  Angle  Corrections 


Lona-Rancre  Case  (R  300  n.  mi.) 

00=  5° 

00=200 

Initial  range  bias  R^ 

(ft) 

75 

22 

Residual  range  bias  <T 

(ft) 

0.75 

0.2 

%  residual  error 

1 

1 

Initial  angle  bias 

(urad) 

3500 

900 

Residual  angle  bias  0 

(urad) 

70 

20 

%  residual  error 

2 

2 

Short-Ranqe  Case  (R  ®  50  n. 

mi.) 

Initial  range  bias  R^ 

(ft) 

22 

7 

Residual  range  biaa<r^jj 

Ut) 

0.5 

0.15 

%  residual  error 

2 

2 

Initial  angle  bias 

(urad) 

2000 

700 

Residual  angle  bias 

(urad) 

60 

20 

%  residual  error 

3 

3 

(Values  shovm  should  be  doubled  for  disturbed  meteorological 
conditions  such  as  heavy  cloud  cover,  fronts,  and  inversions; 
also  for  lack  of  reliable  and  frequent  soundings  covering 
the  entire  tropospheric  path  used  in  measurement) . 
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.  EFFECT  OF  THE  ATMOSPHERE  ON  ACCURACY  OF 

INSTRUMENTATION  SYSTEMS 

4.1  Relationship  between  Position  and  Velocity  Errors. 

Since  the  atmosphere  introduces  errors  in  measurement  of 
target  position,  it  will  also  cause  errors  in  target  velocity 
and  higher  derivatives  of  position.  In  order  to  evaluate  the 
atmospheric  effects  on  velocity  error,  it  has  been  found  con¬ 
venient  to  divide  the  problem  into  two  parts.  The  first, 
which  was  discussed  in  Section  2.  of  this  report,  is  the  effect 
of  temporal  fluctuations  in  the  atmospheric  paths  between  the 
measuring  instruments  and  a  fixed  point  in  space.  The  rate 
of  these  fluctuations  is  governed  by  wind  speed  and  the  presence 
of  turbulence  within  the  atmosphere  through  which  the  measure¬ 
ment  is  made.  This  "direct"  component  of  velocity  error  is 
independent  of  the  target  velocity. 

A  second  term,  which  may  be  described  as  "indirect",  is 
due  to  the  motion  of  the  measurement  paths  through  the 
atmosphere  as  they  follow  a  moving  target,  and  is  dependent 
upon  the  target  velocity  rather  than  on  the  velocity  of  the 
atmosphere  relative  to  the  earth.  The  two  effects  may  be 
combined  in  a  single  term,  but  since  most  of  the  measured 
data  contains  only  the  "direct"  component  it  must  be  mod¬ 
ified  before  being  applied  to  high-velocity  target  tracking. 

The  problem  here  is  similar  to  that  of  relating  bias  (or 
systematic)  errors  of  target  position  to  velocity  errors . 

Brown  (1958)  pointed  out  a  relationship  that  had  been  over¬ 
looked  by  many  in  the  instrumentation  field: 

"It  is  to  be  emphasized  that,  contrary  to  popular 
opinion,  the  determination  of  extremely  precise  ve¬ 
locities  requires  a  correspondingly  high  degree  of 
absolute  precision  of  observations  of  the  tracking 
system.  This  arises  from  the  fact  that  tracking 
systems  measure  space  coordinates  indirectly  only. 

Consequently  small,  constant  biases  in  the  basic 
observations  are  transformed  into  small,  but  slowly 
varying  biases  in  the  space  coordinates .  This  pro¬ 
duces  an  effect  equivalent  to  a  warpage  or  a  dis¬ 
tortion  in  the  scale  of  the  space  coordinate  system. 

If  an  accuracy  of,  say,  1  part  in  100,000  is  desired 


for  velocity,  the  distances  between  trajectory  points 
must  be  unbiased  to  1  part  in  100 , 000 .  Hence  the  need 
for  high  absolute  accuracies . " 


The  effect  of  the  "indirect"  errors  for  the  moving  target 
is  most  pronounced  in  systems  using  range  and  range-rate  data, 
because  in  these  systems  the  atmospheric  limits  have  been 
most  closely  approached.  It  was  stated  in  Section  3.  of  this 
report  that  the  angle  corrections  were  more  sensitive  to 
variations  in  the  refractivity  profile  because  they  depend 
upon  the  gradients  of  this  profile,  measured  in  both  vert¬ 
ical  and  horizontal  planes.  In  the  moving-target  case,  the 
range  and  radial  velocity  errors  will  also  depend  upon  these 
gradients,  since  the  measurement  ray  will  move  through  hundreds 
or  thousands  of  feet  in  the  troposphere  during  the  period  of 
measurement . 

In  calculating  the  effect  of  motion  of  the  beeun,  we  may 
use  either  of  two  methods,  which  yield  identical  results. 
Starting  from  the  range  error  spectrum  of  Pig.  2-6,  which 
applied  to  a  fixed  path,  we  may  derive  expressions  for  both 
position  and  velocity  errors  over  the  path: 

2  fmax 

'r  -  /  df  (4-1) 


(2>r)2 


W(f)  df 


(4-2) 


Here,  the  spectral  density  of  the  range  error  is  represented 
by  W(f) ,  in  units  of  length  ^/cps,  and  the  errors  are  those 
observed  over  a  period  tj^=l/(2f]^)  .  Integration  of  the  spectra 
shown  in  Fig.  2-6  indicates  that  the  position  error  is  quite 
sensitive  to  the  lower  limit  of  the  integral,  corresponding 
to  the  period  of  the  observation.  The  value  of  4^  varies  from 
about  0.1  foot  over  15  to  30  minutes,  to  almost  ten  feet  over 
a  period  of  a  year,  for  paths  at  6°  to  7°  elevation  extending 
through  the  entire  atmosphere.  On  the  other  hand,  the  velocity 
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error  is  almost  independent  of  the  period,  and  also  of  the 
smoothing  time  (up  to  several  seconds) .  Most  of  the  velocity 
error  appears  in  the  spectral  region  from  0.003  to  0.01  cps, 
which  contributes  about  0.001  ft/sec  to  the  total  error  for 
the  fixed  path. 

To  modify  the  spectrum  of  Fig.  2-6  for  the  moving  path, 
we  must  consider  the  velocities  and  spatial  distributions 
of  the  atmospheric  variations  causing  the  error.  The  wind 
velocities  encountered  during  the  Maui  tests  averaged  about 
ten  ft/sec  across  the  path,  and  Pig.  4-1  shows  good  agree¬ 
ment  between  the  spectral  data  and  that  representing  the 
spatial  correlation  of  error,  derived  from  Pig.  2-7.  If 
the  spectrum  were  solely  the  result  of  motion  of  the  air 
mass  across  the  path,  a  frequency  of  0.01  cps  would  cor¬ 
respond  to  a  distance  of  1000  feet  at  ten  ft/sec. 

When  the  target  being  tracked  has  a  tangential  velocity 
component  (normal  to  the  measurement  path) ,  the  relative 
velocity  of  the  troposphere  with  respect  to  the  measurement' 
path  is  increased  by  an  amount  which  varies  with  distance 
from  the  instrument.  Using,  as  an  effective  velocity,  the 
beam  velocity  at  an  altitude  equal  to  one-half  the  scale 
height,  or  at  12,000  feet,  the  relative  velocity  v^  due  to 
beam  motion  will  be  (Fig.  4-2) : 


=  T 
2  s 


(i>  =  i.  s 

2  TT 


(4-3) 


Here,  Rg  is  the  range  to  the  point  where  the  beam  reaches 
the  scale  height,  R  is  the  total  range  to  the  target,  and 
v^  is  the  tangential  velocity  of  the  target.  For  example, 
in  a  satellite  track,  with  R=660  miles  (4x10^  ft) ,  elevation 
0q=6°,  and  v^=10,000  ft/sec,  we  would  have: 

Rg=  210,000  feet,  <4^2.5  mr/sec,  Vjj=260  ft/sec 
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vsloeit! 


Angular  vulocity 


f ,  - > 


b.  BUMB  aotlm  In  horlsontal  plana. 


fi9axa.4-2.  Nation  of  Naaauxaaont  BoaN  for  Noving  Sargat 


In  this  case  the  spectral  distribution  plotted  in  Pig.  2-6 
will  be  shifted  upwards  in  frequency  by  a  factor  of  26,  and 
the  velocity  error  0^  will  be  increased  by  the  same  factor, 
to  about  0.026  ft/sec.  The  predominant  spectral  components 
will  now  be  from  0.08  to  0.25  cps,  and  some  reduction  in 
error  will  be  possible  using  smoothing  times  of  a  few  seconds. 

The  second  way  of  describing  the  "indirect"  component  of 
velocity  error  is  to  consider  the  geometry  of  the  measure¬ 
ment  path,  and  the  exact  direction  of  the  path  as  it  passes 
the  target.  The  radial  velocity  observed  at  the  instrument 
(by  Doppler  or  differentiated  range  measurements)  represents 
the  projection  of  the  target  velocity  vector  on  the  measure¬ 
ment  ray  path  at  the  target.  In  traversing  the  atmosphere, 
the  ray  is  bent  through  an  angle  At*  relative  to  the  assumed 
path,  which  has  been  estimated  using  the  best  available 
tropospheric  profiles. 

It  was  shown  by  Millman  (1958)  that  the  effect  of  ray 
banding  on  Doppler  velocity  measurement  was  more  important 
than  the  direct  error  due  to  the  uncertain  propagation 
velocity  at  the  target  or  along  the  path.  He  expressed  this 
component  of  velocity  error  as  the  product  of  the  target's 
tangential  velocity  times  the  ray  error  angle  at  the  target : 


Av  =  v^.  sin  Aa^  (4-4) 

The  ray  error  angle  describes  the  difference  in  angle  of 
arrival  at  the  target  of  the  measurement  ray  relative  to  the 
angle  estimated  by  the  measuring  device.  If  no  correction  is 
applied  for  refraction  effects,  the  error  would  be  largest 
in  the  elevation  plane,  and  would  be  equal  to  the  difference 
between  the  total  amount  of  ray  bending  (designated  7  in 
Pig.  2-18  or  T  in  Fig.  3-1) and  the  elevation  angle  error  at 
the  instrument  (4) . 

When  corrections  are  applied  in  the  measuring  system,  a 
major  part  of  the  tropospheric  portion  of  Aa^  can  be  eliminated 
by  computation  in  the  system.  As  shown  in  Fig.  4-3,  there  will 
remain  a  small  residual  error  due  to  fluctuation  and  bias  error 
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in  both  elevation  and  azimuth  at  the  instrument.  In  the 
presence  of  a  tangential  component  of  target  velocity,  this 
ray  error  will  produce  an  error  in  apparent  radial  velocity 
given  by  4v.  It  should  be  pointed  out  that  the  need  for 
correction  applies  both  to  those  systems  which  use  the 
angular  data  given  by  the  antenna  and  to  those  which  attempt 
to  ignore  these  angles  by  interferometric  and  trilateration 
solutions  for  target  position,  as  long  as  radial  velocity 
measurements  are  used. 

For  the  troposphere,  the  uncorrected  part  of  the  ray  error 
angle  may  be  estimated  using  the  geometrical  relationships 
of  Fig.  4-3b.  Considering  a  tropospheric  scale  height  H=7  km 
(23,000  feet) ,  the  beam  will  extend  for  a  range  Rg=H/sin  6q 
below  this  height,  and  the  angle  dot  will  be  related  to  the 
elevation  error  6  by: 


R 

Act.  =  ^  {R>R  ,  flat  earth)  (4-5) 

R,+R 

Similarly,  if  the  residual  elevation  error  after  correction 
(or  the  unpredictable  azimuth  error)  is  ,  the  unpredict¬ 
able  component  of  will  be ;  ^ 


Rg+R 


(4-6) 


For  example,  let  us  assume  that  the  target  is  at  an  altitude 
of  160  miles  (10®  ft)  with  a  range  of  660  miles  (4xl0®ft) , 
corresponding  to  an  elevation  angle  6^=6^.  The  tropospheric 
path  length  Rs  will  be  about  1/20  of  the  total  path,  and  the 
ray  error  will  be  1/20  of  the  uncorrected  elevation  or  azi¬ 
muth  error.  Placing  at  about  50  ^iradians,  this  would 

leave  a  ray  error  residue  of  about  2.5  |iradians.  For  a 

missile  with  a  tangential  velocity  of  10,000  feet  per  second. 


the  radial  velocity  measxured  by  even  the  most  precise  Doppler 
system  would  be  in  error  by  .025  ft/sec.  (the  saune  as  found 
by  the  first  method) .  Fortunately,  this  error  will  be 
reduced  at  higher  altitudes  by  the  factor  Rg/(R+Rg) ,  but  the 


i: 

possible  increase  in  the  tangential  velocity  component  to 
values  near  30,000  feet  per  second  will  lead  to  probable 
values  of  several  hundredths  of  a  foot  per  second  in  most 
radial  measurements  on  objects  in  orbit.  These  errors  will 
consist  of  both  bias  and  fluctuating  components,  as  mentioned 
above,  and  will  be  largely  uncorrelated  over  distances  of 
several  thousand  feet.  This  factor  will  place  a  severe  limit 
on  the  attainable  accuracy  of  long-baseline  interferometer 
systems  when  operating  against  targets  with  substantial 
tangential  velocity  components.  Although  not  too  serious 
in  some  guidance  applications,  where  the  missile  at  fuel 
cut-off  is  moving  along  a  near-radial  course,  it  will  cause 
problems  in  most  other  tracking  missions . 

The  ionospheric  component  of  ray  error  may  be  calculated 
in  similar  fashion,  except  that  the  bending  occurs  at  a  point 
much  nearer  the  target,  and  correction  is  virtually  impossible. 
Thus,  the  ray  error  will  be  approximately: 

Ac*.  «  N^/tan  e.  (4-7) 

Here,  0.  is  the  ionospheric  elevation  angle,  having  a  minimum 
value  or  about  9°  for  rays  which  leave  the  earth  at  an  eleva¬ 
tion  angle  of  zero.  A  plot  of  vs.  frequency  for  a  target 

altitude  of  300  km  is  shown  in  Fig.  4-4.  Two  ground  elevation 

angles  are  shown,  and  the  top  scale  indicates  radial  velocity 
error  in  feet  per  second  for  a  tangential  component  of  10,000 
feet  per  second.  When  a  lower  limit  of  3000  mcps .  is  observed 
for  range-rate  measureing  systems,  the  velocity  error  will  be 
between  0.1  and  0.2  feet  per  second  for  v^=10,000  feet  per 
second,  rising  to  a  maximum  of  three  times  this  value  for 
orbital  objects  crossing  the  beam  at  30,000  feet  per  second. 
Operation  below  3000  mcps.  is  permissible  only  at  night, 
when  high  tangential  velocities  are  not  expected,  or  when 
the  range -rate  accuracy  need  not  be  better  than  about  one 
foot  per  second.  Under  extreme  ionospheric  conditions,  it 
may  be  necessary  to  use  frequencies  well  above  3000  mcps . 
in  order  to  reduce  this  component  of  error  to  tolerable  levels. 

The  effect  of  the  atmosphere  on  the  operation  of  any 
instrumentation  system  can  be  visualized  as  a  drift  in  the 
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Flgurs  4-4.  tmy  trxor  Aagl*  opixstiag 


position  of  the  measuring  instrument .  After  correction  for 
measured  tropospheric  profiles,  the  "virtual  source"  (see 
Pig.  4-3)  will  be  found  to  drift  over  a  small  ellipsoid,  as 
shown  in  Fig.  4-5a.  The  extent  of  the  uncorrected  fluctua¬ 
tions  in  a  period  of  fifteen  to  thirty  minutes  is  about 
0.1  foot  along  the  path  and  perhaps  ten  feet  normal  to  the 
path.  VHien  two  instruments  are  used  in  an  interferometer 
system,  the  drifts  will  be  partially  correlated,  as  shown  by 
Fig.  4-1  and  4-5b.  The  long-period  terms  can  be  reduced  by 
obtaining  tropospheric  refractivity  profiles  at  hourly 
(or  shorter)  intervals  and  at  several  locations  around  the 
instrumentation  complex.  However,  the  location  of  the  virtual 
sources  to  absolute  accuracies  of  one  foot  or  better  appears 
very  questionable,  as  this  is  equivalent  to  correction  of 
angle-tracker  data  to  within  5.0  ^radians.  Ionospheric  errors 
will  remain  important  at  all  frequencies  balow  5000  or  6000 
mcps,  when  radial  velocity  data  approaching  0.1  ft/sec  is 
required . 
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a.  Drift  in  Location  of  Virtual  Source  of 
Maaaureiaant  Ray 


Differential 


I  I 

I  Baseline  i 
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b.  Correlated  and  Uncorrelated  Portions  of  Drift 
for  Interferometer 

Figure  4-5.  Atmospheric  Error  in  Terms  of  Drift  In 


Location  of  Instruments 


4.2  Sensitivity  of  Present  and  Proposed  Systems  to 
Atmospheric  Errors. _ 


In  evaluating  atmospheric  errors,  we  may  consider 
separately  the  three  major  types  of  precision  instrumenta¬ 
tion  system.  These  are  the  trackers  (e.g.  pulsed  radar), 
the  interferometers  (e.g.  Azusa  or  Mistreun)  ,  and  the 
trilateration  systems  (e.g.  Glotrack,  the  Goddard  "Range, 
Range-Rate  System",  or  the  radar  trilateration  system).  The 
errors  to  be  evaluated  for  each  system  are  shown  in  Table 
4-1.  Using  the  data  presented  in  Section  2.  for  the  tropo¬ 
spheric  and  ionospheric  errors ,  and  the  correction  methods 
described  in  Section  3  for  the  troposphere,  it  is  possible 
to  estimate  the  atmospheric  components  of  error  in  both 
position  and  velocity  of  targets  for  each  type  of  system. 

This  will  be  done  below  for  the  case  of  high-altitude  targets 
such  as  satellites,  and  the  principal  sources  of  atmospheric 
error  will  be  identified. 

First,  however,  there  are  some  general  statements  which 
may  be  made  regarding  the  errors  introduced  in  the  various 
systems . 

a)  .  The  error  in  range  measurement,  due  to  the  atmo¬ 
sphere,  is  independent  of  the  type  of  system  and  dependent 
upon  frequency  and  atmospheric  path  length.  Residual  errors 
in  the  order  of  one  foot  are  to  be  expected  for  systems 
operating  between  5000  and  10,000  mcps. 

b)  .  Trackers  and  short-baseline  interferometers 
(b<1000  ft)  will  face  limits  of  angular  accuracy  near  0.05 
milliradians ,  due  both  to  fluctuations  in  the  troposphere 
and  to  residual  bias  errors.  Operating  frequencies  above 
1000  mcps.  are  necessary  to  approach  this  limit. 

c)  .  Systems  using  longer  baselines  will  make  possible 
the  achievement  of  better  accuracies .  The  angular  accuracy 
of  an  interferometer  is  limited  by  the  ratio  of  the  range- 
difference  error  to  the  effective  baseline  b',  measured 
normal  to  the  direction  of  the  target.  Fluctuating  com¬ 
ponents  of  show  a  gradual  increase  from  about  0.001  foot 
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feet  (see  Fig.  2-7)  . 


for  b=10  feet  to  about  0.1  foot  for  b>10^ 

d)  .  Beyond  the  baseline  length  b=10^  feet,  the  short¬ 

term  fluctuating  component  of  <T^r  remains  near  0.1  foot, 
and  system  precision  improves  directly  with  baseline  length. 
The  trilateration  systems  represent  a  limiting  case;  where 
the  range -difference  error  is  equal  to  times  the  individual 
range  error.  Residual  bias  (long-term)  errors  are  to  be  ex¬ 
pected  in  both  interferometer  and  trilateration  systems , 
approaching  foot  for  large  station  separations.  The 

residual  bias  for  separations  in  the  order  of  10^  feet  has 
not  been  established,  but  presumably  lies  between  the  values 
of  the  fluctuating  component  (0.1  foot)  and  the  bias  for 
large  separations  (1  foot) .  All  the  cited  errors  are  typical 
of  cloudy  weather,  for  paths  at  7  to  10  degrees  elevation 
angle . 

e)  .  The  fluctuating  component  of  range-rate  difference 
(see  Fig.  2-8)  approaches  0.001  foot  per  second  for  b=10^, 
assuming  20-second  smoothing  and  a  fixed  target.  The  low 
value  of  this  error  in  rate  is  due  to  the  long  periods  in¬ 
volved  in  the  tropospheric  anomalies.  For  the  systems  using 
long  baselines  or  trilateration,  a  much  larger  error  can  be 
expected  to  result  from  the  effect  of  ray  motion  in  following 
a  moving  target,  discussed  in  Section  4.1.  Typical  values 
for  this  component  will  be  given  for  particular  target  con¬ 
ditions  in  the  next  section. 
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TABLE  4-1 


Identification  of  Propagation  Error  Components 
for  Instrxunentation  Systems 

1.  Trackers  (e.g.  Radar) 

Range  error 
Angle  error  <*g 

2 .  Interferometers  (e.g.  Azusa,  Mistram) 

Range  error 

Range-difference  error  0  or 

dr 

Equivalent  angle  errors  =  /(b  sin  B)  or 

6  dr 

Direction  cosine  error  a  ~  c  /b 

c  4r 

3.  Long-Baseline  Trilateration  (e.g.  Glotrack,  Goddard 

Range,  Range-Rate  System,  or  Radar  Trilateration) 
Range  error 

Target  position  error  o  -  a  x  GDOPxV^  or 

P  r 

Equivalent  angle  error  (at  long  range) 

-^2  sin  p) 

All  error  components  consist  of  bias  and  fluctuating  parts. 
Evaluation  of  first  derivatives  is  usually  necessary. 


a. 


Tracker  Error  Analysis 


The  procedure  for  analysis  of  errors  in  the  tracker 
systems  (e.  g.  radar)  is  straightforward.  For  a  given  target 
range  and  altitude,  the  elevation  angle  0q  and  the  tropo¬ 
spheric  path  length  Rg  are  first  determined.  The  initial 
values  of  tropospheric  range  and  angle  bias  may  be  estimated 
from  Figs.  2-2  and  2-4.  Assuming  that  a  correction  procedure 
as  described  in  Section  3.  is  to  be  applied  to  the  data,  the 
residual  bias  errors  may  be  estimated  at  one-half  to  one 
percent  of  the  Initial  range  bias  and  two  to  four  percent 
of  the  angle  bias  5 .  Fluctuating  components  may  be  estimated 
from  Fig.  2-9,  using  an  effective  aperture  width  of  the 
antenna  in  place  of  the  separation  baseline,  and  finding  the 
range  fluctuation  for  the  single-path  case  as : 


Limit 

b-^oo 


VT 


The  tropospheric  rate  errors  are  somewhat  more  difficult 
to  evaluate.  Smoothed  values  of  range  rate  fluctuation  may 
be  read  from  the  limiting  value  of  •  in  Fig.  2-8.  To  these 
must  be  added  the  errors  due  to  motion  qf  the  ray  in  the 
atmosphere,  including  both  bias  and  fluctuating  components. 
When  smoothed,  the  fluctuating  effects  will  tend  towards  zero 
velocity  error,  but  the  bias  portion  will  change  very  slowly 
as  the  target  passes  through  a  wide  interval  in  elevation  or 
azimuth  angle  and  will  not  be  smoothed  out.  The  constant 
part  of  this  error  can  be  considered  as  a  rotation  of  the 
apparent  coordinate  system  with  respect  to  the  true  coordinate 
system  in  which  the  data  is  expressed,  the  rotation  having  a 
value  equal  to  the  average  residual  bias  in  the  tracker  ele¬ 
vation  angle.  The  angular  rate  errors  may  be  read  directly 
from  Fig.  2-8,  using  the  effective  aperture.  Unsmoothed 
values  will  be  considerably  higher  than  those  shown  in  this 
figure  (see  Table  2-4) ,  and  the  eunplitude  and  spectral  spread 
of  the  rate  errors  will  be  increased  by  motion  of  the  beam 
through  the  atmosphere  in  tracking  a  moving  target. 

Ionospheric  errors  for  the  tracker  are  found  in  the 
same  way  as  the  tropospheric  errors,  except  that  initial 
bias  values  are  used  without  allowance  for  correction 
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(see  Figs.  2-15  and  2-16).  Fluctuating  and  ray  bending 
components  are  found  from  Figs.  2-25  and  4-2.  In  general, 
the  bending  angle  d  a ^  will  contribute  the  major  part  of 
the  effective  range-rate  bias  error. 

Total  RMS  position  and  velocity  errors  for  the  tracker 
are  found  by  summing  (in  an  rms  fashion)  the  several  com¬ 
ponents  of  range  and  angle  errors,  multiplying  the  angle 
errors  by  target  range,  and  summing  the  three  final  compon¬ 
ents  corresponding  to  range,  azimuth  and  elevation  coordin¬ 
ates.  (See  Table  4-2). 

b.  Interferometer  Error  Analysis. 

The  procedure  used  in  analysis  of  interferometer  errors 
is  somewhat  more  complex  than  for  the  tracker.  In  the  first 
place,  two  different  values  of  baseline  are  used,  one  for 
azimuth  and  one  for  elevation  (assuming  elevation  angles  below 
90°).  For  distant  targets,  the  average  elevation  angle  ^  at 
which  rays  arrive  at  the  two  stations  may  be  used  to  find  an 
effective  baseline  b '  for  elevation  measurement : 

b'  =  b  sin 

In  most  cases,  tracking  will  be  carried  out  at  elevation 
angles  between  five  and  twenty  degrees,  with  the  effective 
elevation  baseline  shrinking  to  one-third  to  one-twelfth 
the  azimuth  baseline.  As  a  result,  the  azimuth  error  will 
contribute  very  little  to  the  final  RMS  target  position  and 
velocity  error  compared  to  the  elevation  errors .  Using  the 
effective  elevation  baseline.  Figs.  2-7  and  2-8  will  yield 
fluctuating  components  of  position  and  smoothed  velocity  for 
the  target  with  low  tangential  velocity. 

An  angle  bias  error  must  be  found,  either  by  assuming 
corrected  values  of  the  initial  bias  shown  in  Fig.  2-5  or 
by  estimating  the  residual  bias  in  the  range-difference 
measurement.  An  upper  limit  for  interferometer  angle  bias 
will  be  the  value  of  residual  bias  assumed  for  the  tracker. 
This  corresponds  to  correction  of  the  interferometer  curve 
of  Fig.  2-5  to  leave  a  residue  of  about  20%  of  the  initial 
error.  A  lower  limit  would  be  to  assume  a  residual  bias 
error  in  range-difference  equal  to  the  fluctuating  component. 
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The  former  method  will  be  more  accurate  for  the  short-base- 
line  systems,  while  the  latter  will  be  assumed  here  for 
the  long-baseline  systems. 

The  evaluation  of  velocity  errors  is  made  more  complex 
by  the  interaction  between  error  spectra  and  smoothing  time 
for  the  system.  Instantaneous  rate  errors  will  be  governed 
primarily  by  the  product  of  the  ray-bending  error  (see 
Eq.  4-6)  and  the  tangential  velocity.  When  smoothing  is 
employed,  the  errors  will  be  reduced,  as  shown  in  Fig.  4-6. 
Calculations  were  made  from  the  spectrum  of  Fig.  2-6 
(v];jslO  ft/sec)  and  from  a  similar  spectrum  shifted  upwards 
in  frequency  by  a  factor  of  twenty-five  {vjj=250  ft/sec)  . 

The  spectral  density  was  adjusted  to  obtain  the  same  total 
range  error ,  but  the  shift  in  the  spectral  location  of  the 
errors  was  found  to  cause  much  larger  velocitv  errors,  as 
would  be  expected.  Two  curves  are  shown  for  the  high- 
velocity  target  case,  one  with  the  20-second  smoothing  used 
in  Fig.  2-8,  and  one  with  a  2-second  smoothing  time,  as 
often  used  in  range  instrumentation.  The  latter  curve 
shows  the  expected  increase  of  approximately  25:1  over  the 
curve  for  low  velocity.  It  should  be  noted  that  the  curve 
for  2 -second  smoothing  levels  out  sharply  below  b*320  feet 
when  large  apertures  are  used,  due  to  the  smoothing  effect 
of  the  aperture  on  the  received  ray.  When  a  thirty-foot 
aperture  is  used,  the  error  is  35  urad/sec.  The  curves 
plotted  in  Fig.  4-6  are  for  maximum  errors,  and  values  about 
70%  as  great  will  be  used  for  the  "average  weather"  case  to 
be  discussed. 


c .  Trilateration  Error  Analysis 

The  analysis  of  errors  in  the  trilateration  system  is 
quite  simple,  since  all  range  errors  may  be  assumed  equal 
and  uncorrelated.  An  equivalent  angle  error  will  be  found 
according  to  Table  4-1,  assuming  a  GDOP  factor  of  three  for 
the  system.  This  is  equivalent  to  assuming  that  the  measure¬ 
ment  rays  intersect  at  an  angle  near  20°  or  160°,  as  would 
be  the  case  if  a  satellite  were  being  tracked  from  two 
stations  separated  by  about  1000  miles,  with  the  satellite 
midway  between  at  low  altitude.  The  analysis  of  velocity 
errors  will  be  based  on  the  ray  bending  theory,  with  no 
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allowance  for  smoothing.  As  shown  in  Fig.  4-6,  the  effect 
of  smoothing  is  relatively  slight  for  very  long  baselines, 
since  the  significant  errors  have  very  long  periods. 
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4.3  Results  for  High-Altitude  Targets . 


(  ) 


The  foregoing  data  and  analysis  will  be  applied  to  the 
case  of  a  satellite  in  orbit,  at  an  altitude  of  160  nautical 
miles  (10®  feet) ,  to  illustrate  the  problems  imposed  by  the 
atmosphere  on  precision  tracking  systems.  Three  different 
configurations  will  be  assumed,  with  two  operating  frequencies 
compared.  The  results  are  siimmarized  in  Tables  4>2  through 
4-4.  In  each  case,  the  target  parameters  are  as  follows: 

Target  altitude:  h  =  160  n.mi.  =  10^  ft 

Target  range:  R  =  660  n.mi.  =  4xl0°ft 

Tangential  velocity:  Vt=  10,000  ft/sec 

Angular  rate  of  beam  motion :  u)  -  2 .5  mr/sec 

Effective  tropospheric  beam  velocity:  V'jj=  250  ft/sec 
Beam  elevation  angle:  ©o  =  6°  =  105  mr 

The  errors  will  be  given  for  "average"  weather  conditions, 
corresponding  to  the  median  curves  of  Pigs.  2-7  and  2-8,  or 
to  presence  of  small  amounts  of  cumulus  clouds .  The  errors 
for  heavy  cloud  coverage  will  be  about  twice  the  values 
shown,  while  for  clear  sky  they  will  be  about  one-half  as 
great.  Similarly,  the  ionosphere  will  be  assumed  to  follow 
the  characteristics  of  the  daytime  model  (Fig.  2-14) ,  without 
extreme  sunspot  conditions  or  other  "disturbed"  characteristics. 
Where  the  ionospheric  errors  are  important,  the  variation  in 
level  may  be  taken  as  about  three  to  one,  both  above  and  below 
the  values  given,  depending  upon  time  of  day  and  sunspot  cycle. 

The  velocity  errors  for  the  tracker  and  interferometer 
systems  have  been  found  from  Figs.  2-7,  2-8  and  4-6,  taking 
into  account  the  measured  spectra  of  tropospheric  errors. 

For  the  trilateration  system,  the  predominant  error  has  been 
calculated  from  the  uncertainty  of  the  ray-bending  component . 
Total  error  has  been  found  as  the  rms  sum  of  all  bias  and 
fluctuating  components,  and  has  been  expressed  in  terms  of 
range,  angle,  target  position  and  target  velocity.  In  each 
case,  the  components  due  to  angular  measurement  (or  equivalent) 
are  seen  to  govern  the  system  accuracy,  and  of  these,  the 
elevation  component  is  of  greatest  importance .  The  results 
are  not  consistent  with  some  of  the  published  figures  for 
interferometer  systems,  but  are  believed  to  represent  the 
most  accurate  values  for  this  tracking  problem.  The  primary 
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cause  of  the  difference  lies  in  the  fact  that  the  satellite 
may  have  appreciable  tangential  velocity,  causing  the  measure¬ 
ment  beams  to  move  through  the  atmosphere  at  rates  which 
greatly  magnify  the  frequencies  in  the  atmospheric  error 
spectrum.  If  similar  calculations  were  made  for  targets 
which  had  little  or  no  tangential  velocity  (  as,  for  instance, 
missiles  traveling  directly  away  from  the  instrumentation  site) 
the  errors  would  be  appreciably  lower  in  magnitude. 

Table  4-2 .  Typical  Atmospheric  Errors 
Tracking  Radar  on  Satellite  Tracking  Mission 


A.  Tropospheric  Components  (average  weather) 


Range  bias  dR^  (Fig.  2-2)  : 

60  ft 

Angle  Bias  A  (Fig.  2-4)  : 

2500  fir  ad 

Residual  range  bias  (Table  3-1)  : 

0.3  ft 

Residual  angle  bias  (Table  3-1)  : 

50  firad 

Range  fluctuation  <fj.f  (Fig.  2-9)  : 

0.1  ft 

Angle  fluctuation  OQf  (Fig.  2-7)  : 

60  firad 

Range  rate  bias  <T.  =v  <1,  : 

rb  t  Bb 

0.5  ft/sec 

Range  rate  fluctuation  0.  =v  <7  : 

rf  t  ef 

0.6  ft/sec 

Angle  rate  fluctuation  0.  (Fig.  4-6)  : 

0 

4  firad/sec 

B.  Ionospheric  Components  (normal  ionosphere) 


Operating  frequency 

2000 

6000  mcos 

Range  bias  dr^  (Fig.  2-15) 

: 

10 

1.1  ft 

Range  fluctuation  0^.  (Fig. 
Angle  bias  6^  ^(Fig. 

2-25)  : 

0.05 

.006  ft 

2-16) : 

8 

0.9  firad 

Ray  error  4a ^  (Fig. 

4-4)  ; 

30 

3 . 3  firad 

Range  rate  error  0.  =v.da 

r .  t  j. 

0.6 

.066  ft/sec 

i 
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C.  Total  Error 


Operating  frequency 

2000 

6000  mcps 

Range  error  <T^ 

10 

1.2  ft 

Angle  error  0^ 

78 

78  flrad 

RMS  target  position 

310 

310  ft 

RMS  target  velocity  <r^ 

16 

16  ft/sec 

Significant  error  components : 

Residual  tropospheric  range  and  angle  bias 
Ionospheric  range  bias  (2000  mcps  only) 


Table  4-3 .  Typical  Atmospheric  Errors 
Interferometer  System  (Mistram)  on  Satellite  Track 

A.  Tropospheric  Components  (average  weather) 


Range  bias  ^ Rg  (Pig.  2- 

-2)  ! 

60  ft 

Angle  bias  Afi'  (Fig.  2- 

-5)  ; 

180  (irad 

Range -difference  bias  (El.  b'=1000  ft) 

0.18  ft 

Residual  range  bias 

(Table  3-1)  ; 

0.3  ft 

Range  fluctuation 

(Fig.  2-9)  : 

0.1  ft 

Elevation 

Azimuth 

Effective  baseline  b' 

(position  data) : 

1000 

10,000  ft 

Residual  range -difference 

.02 

.06  ft 

Residual  angle  bias  ^"0^: 

20 

6  /irad 

Range -difference 

fluctuation 

2-7)  .014 

.C6  ft 

Angle  fluctuation 

(Pig.  2-7) : 

14 

6  ^rad 
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Table  4-3  (continued) 


Effective  baseline  b* 
(velocity  data)  : 

Elevation 

10,000 

Azimuth 

100,000  ft 

Range -rate  difference 

fluctuation  (T. .  : 

art 

.006 

.0014  ft/sec 

Angle  rate  fluctuation 

(Pig-  4-6)  : 
or 

0.6 

0.14  (irad/sec 

B. 

lonosoheric  Comoonents  (normal  ionosohere.  f=10.000  me) 

Range  bias  ^ri  (Pig.  2-15) : 

0.2 

ft 

Range  fluctuation  (Pig. 

2-15)  :  .001 

ft 

Angle  bias  (Pig.  2-16) : 

0.16  firad 

Ray  error  (Pig.  4-4)  : 

1.0 

firad 

Range  rate  error  : 

.01 

ft/sec 

Ray  difference  error 

0.1 

0.3  firad 

Range  rate  difference  error 

.001 

.003  ft/sec 

Angle  rate  error 

0.1 

0.03  firad/aec 

C. 

Total  Error 

Range  error 

0.35  ft 

Angle  error  <r^ 

25  firad 

RMS  target  position 

100  ft 

RMS  target  velocity 

2 . 4  ft/sec 

Significant  error  sources : 
Tropospheric  range  bias 
Tropospheric  angle  bias  and 

fluctuation 
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Table  4-4.  Typical  Atmospheric  Errors 
Wlde-Basellne  Trllateratlon  System  on  Satellite  Track 


)spherlc  Components 


(average  weather) 


Range  bias  AR^  (Fig.  2-2) : 

Residual  range  bias  (Table  3-1)  : 

Range  fluctuation  (Fig.  2-9)  : 

Geometrical  dilution  factor: 
Equivalent  angle  bias 
Equivalent  angle  fluctuation 
Ray  error  bias 
Range  rate  bias 

Equivalent  angle  rate  bias  0*^^: 

Ray  error  fluctuation 
Range-rate  fluctuation 

Equivalent  angle  rate  fluctuation 


60  ft 
0.3  ft 
0.1  ft 
3 

0.33  ittrad 
0.11  ^rad 
2.5  /irad 
.025  ft /sec 
.027  rad/sec 

1.4  Hr&d 
.014  ft/sec 

.015  /irad/sec 


Ionospheric  Components 


(normal  Ionosphere) 


2000  mcps 

6000  mcps 

Range  bias  Ar^  (Fig.  2-15)  : 

4.5 

0.5  ft 

Range  fluctuation 
(Fig.  2-25)  : 

.03 

.0033  ft 

Equivalent  angle  bias 

4.8 

0.53  /irad 

Equivalent  angle  fluctuation 

0f  • 

.032 

.0035  ^rad 

Ray  error  bias  ^  (Fig .  4-4)  : 

20 

2.2  11  rad 

Range-rate  bias 

0.2 

.022  ft /sec 

Equivalent  angle  rate  bias 

0.22 

.024  ^rad/sec 
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Table  4-4  (continued) 


C.  Total  Error 
Frequency 

Range  error  <r^ 

Equivalent  angle  error  <^0 
RMS  target  position 
RMS  target  velocity  O'y 


2000  mcps  6000  mcps 


4.5  0.6  ft 

4.8  0.64  /trad 

19  2 . 5  f t 

0.9  0.15  ft/sec 


Significant  error  components; 


At  2000  mcps : 
At  6000  mcps : 


Ionospheric  range  bias  and  ray  error 
Ionospheric  and  tropospheric  range  bias  and 
ray  error  (tropospheric  is  slightly 
greater  than  ionospheric  at  this  frequency) 
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CONCLUSIONS 


5.1 

Tropospheric  bias  errors  are  highly  predictable  using 
radiosonde  or  refractometer  profiles ;  residual  errors  from 
1?6  to  3%  of  the  initial  bias  levels  are  commonly  attained 
using  procedures  described  in  Section  3 .  Data  to  within  one- 
half  foot  in  range  and  20  to  70  ^radians  in  angle  can  be  ex¬ 
pected  at  elevation  angles  above  five  degrees . 


5.2 


Tropospheric  fluctuation  errors  are  not  correctable 
using  any  known  procedure,  and  will  amount  to  a  few  tenths 
of  a  foot  in  range,  and  10  to  50  /tradians  in  angle  (depend¬ 
ing  on  the  baseline  or  aperture  used  for  measurement) ,  under 
normal  weather  conditions . 


5.3 


The  relationship  between  temporal  and  spatial  correlation 
of  tropospheric  fluctuations  has  been  investigated,  based  on 
data  obtained  by  the  National  Bureau  of  Standards .  The  effect 
of  short-period  fluctuations  is  described  by  Figs.  2-7  and 
2-8,  and  is  consistent  with  a  drift  of  tropospheric  anomalies 
across  a  fixed  measurement  path  at  the  speed  of  the  prevailing 
wind. 


5.4 


In  range  instrumentation  applications,  where  the  beam  is 
not  fixed,  the  residual  "bias"  and  long-term  error  components 
will  change  as  the  beam  moves,  and  additional  atmospheric 
rate  errors  will  be  generated,  as  shown  in  Fig.  4-6.  These 
errors  will  be  proportional  to  the  tangential  velocity  of 
the  missile,  and  will  typically  be  five  to  fifty  times  the 
errors  measured  for  a  fixed  beam. 


5.5 


The  uncertainty  in  tropospheric  path  leads  to  errors 
equivalent  to  motion  of  the  instrument  on  the  ground.  The 


motion  of  the  "virtual  source"  typically  amounts  to  several 
feet  normal  to  the  path  and  a  few  tenths  of  a  foot  along  the 
path. 


5.6 


Ionospheric  errors  are  essentially  unpredictable,  and  will 
exceed  the  residual  tropospheric  errors  when  operating  fre¬ 
quencies  below  3000  mcps  are  used.  Even  in  the  5000-6000  mcps 
band  the  ionospheric  errors  will  contribute  to  overall  atmos¬ 
pheric  error  during  daytime  operation. 


5.7 


Redundant  measurements  performed  at  two  frequencies  below 
3000  mcps  can  be  used  to  correct  for  ionospheric  error  in 
both  range  and  angle. 


5.8 


The  lowest  atmospheric  errors  are  found  in  trilateration 
systems  using  very  long  baselines.  Total  position  and  velocity 
errors  for  a  typical  satellite  track  (660  miles  range,  100 
miles  altitude)  through  average  weather,  are  as  follows: 


RMS  Position 
Error  (feet) 

RMS  Velocity 
Error  (ft/sec) 

Range-angle  tracker  at 
6000  mcps 

310 

16 

Interferometer , 

10,000  mcps  (Mistram) 

100 

2.4 

Trilateration  system, 
2000  mcps 

6000  mcps 

19 

2.5 

0.9 

0.15 

The  above  errors  may  be  increased  or  decreased  by  a 
factor  of  two  or  three  for  different  weather  conditions 
(and  at  2000  mcps  for  different  ionospheric  conditions) . 
The  trilateration  errors  shown  are  dependent  upon  perfect 

survey  of  station  location,  as  well  as  instrumental  errors 
below  one-half  foot  in  range  and  0.02  ft/sec  in  range-rate. 


6.  RECOMMENDATIONS  FOR  ACHIEVING  INCREASED  ACCURACY  TODAY 


6.1 


Since  ionospheric  refraction  cannot  at  present  be  predic¬ 
ted  to  within  better  than  about  50%  of  any  instantaneous 
value,  the  use  of  microwave  bands  or  of  dual-frequency 
measurements  is  necessary  in  precision  tracking  of  targets 
above  100  miles.  Single  frequency  systems  requiring  velocity 
data  better  than  one  foot  per  second  should  operate  above 
3000  megacycles  to  minimize  the  ionospheric  refraction  effect. 


6.2 


A  continuing  program  of  data  analysis  at  the  various 
ranges  should  be  instituted  to  evaluate  and  improve  the 
atmospheric  correction  procedures  described  in  Section  3.  of 
this  report.  Some  of  this  work  is  being  done  at  AMR  now  (and 
perhaps  at  other  Ranges)  in  connection  with  other  activities . 
However  this  work  is  so  important  that  it  should  be  support¬ 
ed  as  a  separate  function;  this  is  the  only  way  it  will 
receive  the  attention  which  it  deserves .  Data  is  available 
at  all  of  the  Ranges ;  it  is  only  necessary  that  qualified 
people  be  assigned  to  an  analysis  of  it.  This  work  should 
be  fully  supported  by  the  Air  Force  at  its  Ranges  and  by 
the  other  services  at  their  respective  Ranges .  This  data 
analysis  conducted  at  each  particular  Range  should  be  fully 
coordinated  cunong  the  Ranges  and  with  the  NBS  measurement 
program.  Other  interested  agencies,  such  as  NASA,  should 
also  be  invited  to  participate.  The  Inter  Range  Instr\jmen- 
tation  Group  has  done  an  excellent  job  in  the  past  of 
providing  coordination  and  dissemination  of  technical  infor¬ 
mation  among  the  Ranges  on  an  informal  basis  and  is  well 
qualified  to  do  so  in  this  case.  This  coordination  of  effort, 
especially  among  the  National  Ranges ,  should  receive  the  full 
support  of  DOD. 


6.3 


Standard  procedures  for  atmospheric  correction  of 
tracking  data  should  be  adopted  by  all  of  the  Ranges  and 
Range  users  and  estimates  of  residual  bias  errors  agreed 
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upon  for  each  procedure.  The  methods  discussed  in  Section  3 
of  this  Report  are  suggested  as  a  basis  for  such  Standards 
and  are  consistent  with  efforts  now  underway  by  the  Electro¬ 
magnetic  Propagation  Working  Group  of  the  Inter  Range 
Instrxunentation  Group.  The  EPWG  is  currently  working  on  a 
range  instrumentation  manual  which  it  hopes  will  lead  to 
more  standardization.  The  work  of  the  individual  members  of 
EPWG  on  this  project  should  be  given  the  full  support  of 
each  particular  Range  where  they  are  located  and  the  project 
as  a  whole  should  have  the  complete  support  of  DOD  on  an 
inter -range  basis. 
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1. 


RECOMMENDED  RESEARCH  FOR  FUTURE  INCREASE  IN  ACCURACY 


i  > 


7.1 


Future  tracking  systems  should  be  designed  to  tolerate 
the  unpredictable  fluctuations  of  the  measurement  ray  paths 
in  the  atmosphere.  When  targets  of  high  velocity  must  be 
tracked  with  accurate  three-coordinate  velocity  measurements, 
the  measurement  systems  baselines  should  be  as  long  as  possible 
they  should  be  consistent  with  target  altitude.  Systems  which 
require  that  instrxamentation  sites  be  located  with  an  accuracy 
on  the  order  of  one  foot  or  better  do  not  appear  to  be  consist¬ 
ent  with  our  ability  to  predict  the  ray  paths  in  the  tropos¬ 
phere  . 


7.2 


A  specific  procedure  for  measuring  and  correcting  tropo¬ 
spheric  errors  on  a  real-time  basis  has  been  proposed  to  the 
Panel.  A  brief  discussion  of  this  technique  is  given  in 
Appendix  B.  Theoretically  this  technique  appears  very  promis¬ 
ing.  It  is  now  a  question  of  determining  whether  experimental 
verification  can  be  obtained.  This  work  should  receive  full 
support  from  the  Air  Force. 

7.3 


The  National  Bureau  of  Standards  (Boulder)  has  outlined 
a  program  of  atmospheric  measurements  which  it  is  attempting 
or  would  like  to  attempt.  This  prograun  is  discussed  in 
Appendix  C.  These  measurements  would  provide  much  needed 
data  on  spatial  and  temporal  correlation  of  tropospheric 
range  (or  phase)  errors .  The  Panel  recommends  that  this 
prograun  be  pressed  as  rapidly  as  possible  and  fully  supported 
by  the  Air  Force  to  provide  much  needed  information  for  both 
the  interferometer  systems  and  the  longer  baseline  systems 
using  range  and  range -rate  data. 


7.4 


When  tracking  at  interplanetary  distances  the  errors 
imposed  by  the  atmosphere  become  proportionally  less .  Conse¬ 
quently  one  limiting  factor  to  tracking  accuracy  at  such 
distances  would  appear  to  be  the  precision  with  which  we 
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know  the  velocity  of  light,  currently  felt  to  be  about  1  part 
per  million.  (A  discussion  of  our  knowledge  of  the  velocity 
of  light  is  given  in  the  Report  of  the  Ad  Hoc  Panel  on  Basic 
Measurements.)  The  efforts  by  the  National  Bureau  of  Stand¬ 
ards  to  better  determine  this  value  should  be  fully  supported. 


7.5 


With  our  current  tracking  techniques  for  interplanetary 
distances  an  even  more  critical  need  than  a  better  determina¬ 
tion  of  "c"  is  that  of  a  better  frequency  standard.  For 
doppler  tracking  a  target  over  such  distances  a  frequency 
st^^dard  or  clock  having  a  short  time  stability  of  1  part  in 
10‘‘‘  is  needed  now.  (At  the  present  time  we  can  measure  time 
with  an  accuracy  of  about  1  part  in  10^^.  This  is  discussed 
in  the  Report  of  the  Ad  Hoc  Panel  on  Basic  Measurements.)  The 
continuing  efforts  of  the  National  Bureau  of  Standards  to 
develop  more  stable  frequency  standards  should  receive  full 
support . 


7.6 


Since  the  accuracy  limits  of  current  tracking  instruments 
and  propagation  correction  procedures  are  on  the  order  of 
about  one  foot,  there  is  a  definite  need  for  geodetic  systems 
or  procedures  capable  of  locating  our  tracking  instruments 
or  systems  to  this  same  accuracy  over  intercontinental  distances 
It  is  recommended  that  a  study  group  be  convened  to  determine 
what  are  the  most  fruitful  areas  for  investigation  which  could 
lead  to  a  better  determination  of  locations  on  a  global  basis 
and  which  ultimately  perhaps  might  lead  to  the  accuracy  men¬ 
tioned  above. 


{ 
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8 .  FURTHER  RECOMMENDATIONS 


8.1 


The  Ranges  should  make  the  systems  designers  and/or  Range 
users  feuniliar  with  the  basic  limitations  on  tracking 
accuracy  imposed  by  the  atmosphere  as  described  in  this  Report. 
It  is  futile  for  Range  users  to  request  accuracies  which 
cannot  be  obtained  for  reasons  discussed  herein.  And  the 
Panel  does  not  anticipate  that  significant  improvement  over 
the  potential  accuracies  discussed  here  will  be  attained  in 
the  near  future,  although  more  consistent  use  of  correction 
techniques  will  improve  on  past  performance.  However,  if 
the  research  recommended  in  this  Report  is  undertaken  and 
adequately  supported  it  may  disclose  means  of  reducing  the 
basic  uncertainties  connected  with  propagation  through  the 
atmosphere  which  could  be  applied  within  the  next  decade. 


8.2 


In  order  to  make  the  best  use  of  the  available  resources 
for  the  development  of  instrumentation  systems  and  techni¬ 
ques  for  propagation  correction,  the  responsible  agencies 
should  arrive  at  consolidated  requirements  for  missile  and 
satellite  measurement  accuracies  instead  of  new  and  differ¬ 
ent  requirements  for  each  individual  program.  The  consolida¬ 
ted  requirements  should  be  stated  and  published  in  such  a 
manner  as  to  encourage  scientific  work  on  the  most  fundeunental 
instrumentation  problems,  and  should  not  be  hampered  by  the 
security  restrictions  and  "need-to-know"  of  any  particular 
weapon  or  weapon  system.  This  is  a  DOD  wide  problem  and  DOD 
should  take  the  lead  in  trying  to  implement  this .  However , 
the  Air  Force  could  do  much  along  this  line  with  those 
progreuns  under  its  cognizance. 
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Appendix  A 


OBSERVED  ELEVATION  ANGLE  ERRORS 
by 

Dr.  John  B.  Smyth 
Smyth  Research  Associates 


The  smooth  gross  structure  of  the  N  profile  causes  errors 
which  can  be  corrected  if  one  knows  the  shape  of  the  N  profile, 
but  the  inhomogeneities  cause  an  indeterminancy  in  the  actual 
vehicle  location  which  cannot  be  resolved.  The  principle  con¬ 
cern  in  precision  tracking  is  to  correct  for  the  gross  effects 
(the  bias  errors)  and  to  estimate  as  well  as  possible  the 
magnitude  of  the  random  errors. 

Smyth  Research  Associates,  under  contract  to  the  Air  Force*, 
are  using  a  sea  interferometer  technique  to  measure  the  apparent 
position  of  satellites  as  they  arise  above  the  horizon  over  the 
Pacific  Ocean.  This  experimental ' facility  offers  a  high  degree 
of  resolution  for  sources  near  the  horizon.  Figure  A-1  shows 
the  location,  a  picture  of  one  of  the  45'  by  70'  antennas  and 
a  schematic  of  the  experimental  arrangement.  Data  obtained  at 
this  site  have  been  used  to  evaluate  the  refraction  of  the 
ionosphere  and  troposphere  as  a  function  of  elevation  angle. 

Figures  A-2  and  A-3  give  information  on  the  average  error 
angle  at  frequencies  from  150  to  400  Me  for  elevation  angles 
between  zero  emd  ten  degrees.  The  spread  in  the  data  is  several 
milliradlans ,  more-or-less  independent  of  elevation  angle.  The 
reason  for  this  is  not  completely  clear,  although  it  appears 
that  a  large  part  of  this  spread  is  in  the  uncertainty  of  the 
location  of  the  satellite. 

A  comparison  of  observed  and  predicted  elevation  angle 
errors  is  given  in  Figure  A-4. 


*AF30 (602) 2084,  Rome  Air  Development  Center. 


Figure  A-1  Location  of  the  SUA  Point  Buchon  Field  Station. 
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Figure  A-2  Average  error  angle  July  1960. 
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Figure  A-3  Average  error  angle  September  and  November  1961. 


PREDICTED  IONOSPHERIC  ERROR  ANGLE  -  MILLIRADIANS 


Figure  A-4 
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Appendix  B 


A  LINE  INTEGRAL  REPRACTOMETER 
by 

John  F.  Sullivan 
The  Mitre  Corporation 


Introduction 


This  paper  was  prepared  for  the  National  Academy  of  Sciences 
Ad  Hoc  Panel  on  Electromagnetic  Propagation,  and  is  a  report  on 
the  status  of  certain  studies  by  the  MITRE  Corporation  of  the 
ultimate  limitations  on  the  accuracy  and  precision  of  microwave 
measurements  imposed  by  the  troposphere,  and  of  techniques  for 
correcting  for  its  effect.  The  report  is  somewhat  premature 
in  that  no  experimental  confirmation  of  the  analyses  is  avail¬ 
able,  but  since  the  subject  is  the  committee's  "raison  d'etre" 
and  the  conclusions  are  somewhat  at  variance  with  current  prac¬ 
tice  in  this  field,  it  is  felt  that  a  report  should  be  presented 
for  consideration. 

Discussion  of  Problem 


Essentially,  the  conclusion  reached  is  that  a  measurement 
can  be  made  in  real  time,  over  the  path  of  interest,  which  will 
allow  for  the  correction  of  the  effects  of  a  real  turbulent 
atmosphere  on  position  determinations.  Residual  errors  remain¬ 
ing  after  such  a  correction  have  been  estimated  to  be  at  least 
an  order  of  magnitude  down  from  present  limitations  in  accuracy 
due  to  imperfectly  known  inhoroogeneities  and  turbulence. 

To  illustrate  this  consider  a  ray  traversing  a  path  S 
between  points  P^  and  P^  in  a  turbulent,  inhomogeneous  atmos¬ 
phere,  whose  index  of  refraction  is  described  by  some  scalar 
function,  n(x,y,z,t),  of  space  and  time.  The  elemental  time 
(dt)  required  for  a  wave  to  traverse  an  elemental  arc  length 
(ds)  along  the  ray  is  given  by  the  expression: 
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where  C  is  the  velocity  of  light  in  vacuvun. 

It  is  convenient  to  work  with  refractivity  N  defined  by: 

N(s,y,z,t)  =  n(x,y,z,t)  -1  (B-2) 

By  combining  these  expressions  and  integrating  their  results 

S  =  CAt  -  J  N(x,y,z,t)ds  (B.3) 

^  s 

where  At  is  the  transit  time  as  conventionally  measured  by 
electronic  techniques. 


This  expression  indicates  that  in  order  to  know  the  arc  length 
between  two  points  in  such  an  atmosphere , it  is  necessary  to 
measure  the  transit  time,  multiply  this  by  the  vacuum  velocity 
of  light,  and  diminish  the  result  by  a  quantity  which  repre¬ 
sents  the  integrated  effect  of  the  slowing  of  the  wave  by  the 
material  atmosphere. 

It  is  this  latter  quantity  j  N(x,y,z,t)ds  which  we  propose 
to  measure.  Before  discussing  t^e  method  of  measurement,  how¬ 
ever,  it  would  be  useful  to  excunine  the  following  implications 
of  such  a  measurement: 


1.  To  correct  for  the  effect  of  a  turbulent  inhomogeneous 
atmosphere  on  an  estimate  of  arc  length  (ds)  from  a  conventional 
transit  time  measurement  (At)  it  is  only  necessary  to  know  this 
one  time  varying  quantity  Nds) .  Present  correction  tech¬ 
niques  assume  some  simple  average  spatial  functional  relation¬ 
ship  of  N{x,y,z).  More  or  less  sophisticated  attempts  are  made 
to  account  for  inhomogeneities  and  temporal  variations  of  the 
medium  by  using  this  relationship  to  extrapolate  from  one  or 
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interpolate  between  a  relatively  few  measurements  of  point  index 
of  refraction.  The  estimate  of  N{x,y,z,t)  is  based  on  previously 
measured  data  which  becomes  increasingly  more  expensive  to  obtain 
as  the  time  and  space  intervals  between  data  points  is  reduced. 
Prom  this  estimate,  and  its  estimated  gradient,  rays  are  traced, 
and  the  quantity  N(x,y,z,t)ds  is  computed  by  some  more  or 
less  complicated  pfograun.  Here  again,  as  the  time  interval 
between  computer  points  decreases,  the  time  and  expense  of  the 
calculations  required  rapidly  increases.  It  is,  therefore,  at 
least  intuitively  preferable  to  make  a  single  measurement  of 
the  quantity  i  Nds  in  real  time  which  can  be  simply  subtracted, 
again  in  real*'€ime,  from  the  transit  time  measurement  multiplied 
by  the  velocity  of  light.  No  more  sophisticated  computation  is 
necessary. 

Until  some  experimental  data  upon  which  to  base  an  error 
analysis  is  available,  the  above  statements  cannot  be  proven 
but  the  intuitive  argument  at  least  provides  an  incentive  for 
looking  for  such  a  measurement  technique.  We  are  thinking  in 
terms  of  an  exceedingly  accurate  (about  1  part  in  10®)  measure¬ 
ment  of  arc  length.  The  basic  measurement  of  transit  time  is 
more  than  capable  of  this  kind  of  accuracy.  It  should  be  noted, 
however,  that  the  measurement  of  the  correction  necessary  to 
account  for  refraction  need  not  be  made  to  anything  like  this 
kind  of  percentage  accuracy.  Since  refraction  corrections 
typically  cunount  to  a  few  hundred  parts  per  million,  a  measure¬ 
ment  accuracy  of  slightly  better  than  1%  is  adequate. 

2 .  Gradients  in  the  index  of  refraction  structure  have 
the  effect  of  bending  the  ray  path.  This  bending  causes  the 
arc  length  S  to  differ  from  the  straight  line  distance  R 
between  and  P2.  It  can  be  readily  shown  that  this  difference 
in  length  is  completely  negligible  (  <1x10“')  under  normal  at¬ 
mospheric  conditions  unless  the  ray  becomes  trapped  in  a  duct. 

Thus,  a  direct  measurement  of  the  line  integral  of  refrac- 
tivity  can  give  a  significant  improvement  to  measurements  of 
range,  and  angles  derived  from  range  sum  or  range  difference 
measurements,  without  a  detailed  knowledge  of  the  refractive 
index  structure. 

3.  Prom  exeunination  of  the  mathematical  rules  for  the 
differentiation  of  an  integral  whose  limits  are  functions  of 
the  variable  of  differentiation,  one  may  conclude  that  from 
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a  time  derivative  of  this  measurement  N (x,y ,  z ,  t)  ds  one  may 

obtain  a  real  time  correction  for  all  of  the  refractive  effects 
of  the  troposphere  on  the  magnitude  of  doppler  measurements . 

This  correction  would  include  the  effect  of  turbulent  variations 
of  the  medium,  the  effect  of  apparent  variation  in  path  length 
due  to  motion  through  a  inhomogeneous  medium,  as  well  as  the 
effect  of  the  refractive  index  at  the  moving  body.  An  angle 
(a)  between  the  ray  S  and  the  straight  line  P  results  from 
gradients  in  the  refractive  structure.  Prom  thii  measurement 
alone .  however,  one  obtains  no  information  on  the  effects  of 
refraction  on  direction  of  the  ray  at  the  moving  vehicle.  The 
uncertainty  in  this  angle  is  of  significance  in  precision 
velocity  determinations .  Thus ,  for  velocity  measurements ,  this 
measurement  alone  does  not  completely  account  for  refractive 
effects . 

Prom  the  above  considerations,  it  can  be  concluded  that  a 
device  which  would  measure  in  real  time  the  integrated  effect 
of  the  slowing  of  electromagnetic  waves  by  the  troposphere 
could  be  useful  in  significantly  reducing  the  effects  of  present 
limitations  of  knowledge  of  this  region.  This  does  not  mean 
that  there  are  no  further  limitations  to  the  accuracies  attain¬ 
able,  merely  that  they  are  different,  and  from  our  examination 
appear  smaller.  These  include  the  difference  between  arc 
length  and  range  previously  noted  and  the  deviations  from  simple 
ray  theory  upon  which  the  foregoing  discussion  is  based. 

It  would  be  rather  fruitless  to  describe  the  improvements 
possible  by  means  of  a.  measurement  without  some  conviction 
that  such  could  be  accomplished.  We  are  working  on  a  technique 
which  we  believe  will  ultimately  provide  the  desired  infor¬ 
mation.  This  work  is  in  its  early  stages  and  much  remains  to 
be  proven  by  means  of  equipment  design  as  well  as  laboratory 
and  field  tests  of  the  device.  A  number  of  detailed  calcula¬ 
tions  have  been  performed  which  have  led  us  to  the  conclusion 
that  such  an  effort  would  be  fruitful.  While  they  have  no 
place  in  this  report,  we  would  be  happy  to  discuss  them  with 
those  members  who  are  interested  in  the  details.  What  follows 
is  intended  to  be  descriptive  of  the  technique  proposed. 

Proposed  Technique 

The  radio  index  of  refraction  in  any  elemental  volume  of 
the  atmosphere  is  considered  as  the  sum  of  two  parts;  the  dry. 
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or  optical,  part  arising  from  the  atomic  refractivities  of  the 
gaseous  constituents  (N2/  O2,  A,  etc.)  weighted  in  accordance 
to  their  relative  abundance;  and  the  wet  part,  applicable  only 
to  radio  and  microwaves,  arising  from  the  permanent  electric 
dipole  moment  of  the  uncondensed  water  vapor  molecules  present. 
At  sea  level  and  normal  humidity  the  total  index  is  about  300  N 
units  (N  =  (n-1)  x  10®)  with  about  60  N  units  of  this  total 
representing  the  "wet"  contribution.  The  relative  abundance 
of  the  various  dry  components  is  relatively  constant  throughout 
the  atmosphere,  while  the  water  vapor  density  is  quite  erratic. 


In  the  millimeter  wave  region  molecular  reasonances  of  water 
vapor  and  oxygen  occur  producing  high  absorption  in  the  frequency 
region  associated  with  the  spectral  lines .  Because  of  the  in¬ 
terest  in  microwave  propagation,  and  the  basic  information  on 
the  structure  of  the  molecules  inherent  in  its  parameters,  this 
attenuation  has  been  well  explored  both  theoretically  and  ex¬ 
perimentally.  Associated  with  the  resonances  there  is  also  a 
dispersion  or  variation  of  refractive  index  with  frequency. 

These  dispersions  are  extremely  small  (theoretically  calculable 
to  about  .03  N  units  differential  index  across  the  1.35  cm 
water  vapor  line,  and  1.97  N  units  across  the  60  kmc  oxygen 
line,  at  sea  level  and  normal  temperatures  and  humidity) .  These 
values  of  dispersion  are  theoretical  and  have  never  been  experi¬ 
mentally  measured  precisely.  Measurements  by  Essen  of  total 
refractivity  have  been  made  at  9  kmc  and  24  kmc  and  confirmed 
that  the  difference  is  less  than  .12  N  units. 


Assume  that  the  small  differential  index  of  refraction 
between  two  frequencies  (f^^  and  f2  for  HOH,  f^  and  f^  for  O2) 

on  either  side  of  these  spectral  lines  is  proportional  (at 
every  point  in  the  atmosphere)  to  the  contributions  of  the 
wet  and  dry  atmospheric  constituents  to  the  total  refractive 
index.  This  assumption  is  valid  if  the  frequencies  involved 
are  far  enough  away  from  the  resonances  so  that  collision 
broadening  effects  are  negligible  and  temperature  correction 
techniques  (dependent  upon  the  application)  are  emplyed. 

This  assumption  can  be  expressed  as : 


N  -  N  =  K  N 
1  2  HOH  HOH 


(B-4) 
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where 


=  refractivity  at 

N  =  refractivity  at  f 

=  refractivity  at  f^ 
=  refractivity  at  f^ 


HOH,  ^dry  =  proportionality 
constants  under 
assumption 

N 

HOH  =  Wet  component  of 

total  refractivity 

N 

dry  =  Dry  component  of 

total  refractivity 


Consider,  then,  a  system  consisting  of  a  transmitter  which  may 
be  placed  on  a  moving  vehicle  and  a  remote  ground  receiver 
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The  transmitter  radiates  a  CW  wave  which  at  some  point  in  its 
plumbing  may  be  expressed  as : 

A  =  A  sin  (2»rf  t)  +  A  sin  (27lf^  t)  (B-5) 

1  12  2 

where  fj^  and  f^  are  harmonically  related  and  locked  in  phase. 
At  a  remote  receiver  the  phases  of  the  two  components  of  the 
wave  received  may  b«  written  from  the  wave  equation : 

r  n, (t,s)  ds 

*1  =  -  la 

f  n^(t,s)  ds 
- ) 

n]^  is  very  nearly  equal  at  all  points  to  n2  so  that  both 
rays  traverse  essentially  the  same  path.  In  the  receiver  the 
two  components  are  separated,  one  component  is  multiplied  by 
the  frequency  ratio,  and  the  resultant  phases  compared  in  a 
phase  detector.  A  voltage  ®HOH  will.be  produced. 

"hoh  =  ^  ^  (t.s)  -  (t,s) 
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where  is  the  phase  detector  constant. 

But,  by  the  preceding  assvanption, 

(t,s)  -  (t.s)  =  (t,s)  (B-8) 

•  e  K  K  f  ,  ,  ,  , 

.  .  HOH  =  1  HOH  N  (t,s)  ds  (B-9) 

Js  HOH 


If  a  similar  procedure  is  undertaken  at  frequencies  f^  and  f^ 
about  the  Oxygen  line  eind  the  two  voltages  svuraned  the  result 
will  be : 

e  K  f  N  ,  ,  f  N  ,  , 

TOT  =  K  HOH  \  HOH  (t,s)  ds  +  K  K  dry  (t,s)  ds 

1  Js  2  dry  Js  (B-lO) 

e  f  N 

TOT  =  \  TOT  (t,s)  ds  (Upon  proper  adjustment  of  and  K2) 

(B-11) 

This  voltage  is  then  directly  proportional  to  the  correction  to 
the  transit  time  measurement  required  to  account  for  the  effects 
of  refractivity  in  the  atmosphere  over  that  particular  path  and 
at  that  particular  time  required  by  equation  B-3 

It  should  be  noted  that  the  transmiss^n  is  unidirectional. 
The  measurement  is  made  of  the  quantity  J  N(x,y,z,t) .  A  sep¬ 
arate  two-way  transmission  is  required  fo?  the  basic  transit 
time  measurement.  (The  treatment  here  represents  an  attack  on 
the  problem  of  clear  air  inhomogeneities  and  turbulence.  Con¬ 
densed  water  vapor  in  the  form  of  rain  and  clouds  will  atten¬ 
uate  these  frequencies  and,  more  significantly,  lack  the  gaseous 
in  wave  absorption  lines.  Analysis  of  this  latter  problem  is 
continuing,  but  a  solution  to  the  clear  air  problem  seems 
enough  of  a  step  forward  to  implement  the  necessary  equipment.) 

There  are  a  number  of  equipment  problems  associated  with 
the  proposed  measurement  technique  described,  but  a  conclusion 
has  been  made  that  with  care  in  design,  a  measurement  of  the 
extremely  small  phase  differences  can  be  accomplished.  It 
would  be  well  to  state  here,  however,  that  the  reason  that  the 
measurement  is  only  formidable  and  not  impossible  is  the  fact 
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that  it  is  a  differential  measurement.  The  same  atmospheric 
path  is  traversed  by  the  rays  of  both  f and  f  .  Insofar  as 
possible,  f,  and  will  traverse  the  Scune  paths  of  the 
antennae,  plumbing,  transmitter,  and  receiver  elements  which 
will  be  designed  to  provide  phase  shifts  linear  with  frequency 
(calibration  techniques  are  anticipated  to  be  required) .  Doppler 
effects  (neglecting  relativity)  are  proportional  to  frequency 
and  cancel.  Care  must  be  exercised  in  the  design  of  the  antennae 
to  reduce  multipath  to  negligible  proportions  and  preclude  errors 
due  to  slight  misalignment  of  antennae  boresight. 

Conclusion 


An  instrument  can,  therefore,  be  developed  which  will  direct¬ 
ly  measure  the  integrated  effect  of  the  slowing  of  a  radio  wave 
through  a  real  turbulent  inhomogeneous  atmosphere.  This  can  be 
accomplished  by  the  measurement  of  the  small  difference  in  one¬ 
way  transit  time  between  frequencies  on  either  side  of  the 
microwave  water  vapor  and  oxygen  absorption  lines  by  the  phase 
detector  technique  described.  While  a  single  measurement  gives 
no  indication  of  the  bending  of  the  ray  which  has  occured,  this 
can  be  shown  to  be  a  second  order  effect  in  the  estimation  of 
position. 
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Appendix  C 


PROPOSED  EXPERIMENTAL  STUDIES  OF  ATMOSPHERIC 
LIMITATIONS  ON  RADIO  TRACKING  ACCURACY 

by 

Harris  B.  Janes 

National  Bureau  of  Standards,  Boulder 


The  following  is  an  outline  of  some  experimental  research 
that  should  be  undertaken  to  study  the  limitations  imposed  by 
the  troposphere  on  line-of-sight  radio  ranging  and  tracking 
systems.  This  is  not  intended  to  be  an  exhaustive  list  of 
projects.  In  fact  it  is  heavily  weighted  in  the  direction  of 
those  projects  which  constitute  an  extension  of  work  performed 
by  the  National  Bureau  of  Standards  over  a  period  of  several 
years,  and  hence  for  which  specialized  equipment,  techniques 
and  personnel  are  available. 

The  proposed  experimental  work  should  be  divided  into  two 
categories;  "^basic"  and  "applied".  The  basic  atmospheric 
research  consists  of  studying  the  four-dimensional  statistics 
of  the  radio  refractive  index  in  the  turbulent  troposphere. 

This  will  be  done  first  on  a  relatively  modest  scale  with 
microwave  refractometer  sensors  mounted  on  an  orthogonal  three- 
dimensional  array  on  a  very  flat  mesa  near  Boulder,  Colorado. 

The  purpose  is  to  describe  atmospheric  turbulence  in  terms  of 
the  cross  spectra  between  separated  sensors  as  a  function  of 
separation,  and  the  corresponding  spatial  correlation  functions 
and  wave  number  spectra  for  various  conditions  of  wind  velocity, 
solar  radiation,  etc.  Future  extensions  of  this  work  would  in¬ 
clude  flying  an  array  of  refractometer  sensors  on  a  light  air¬ 
craft  to  study  the  spatial  structure  as  a  function  of  altitude. 

It  is  the  second  or  "applied"  research  category  that  will 
probably  be  of  most  immediate  interest  to  people  concerned  with 
missile  tracking  systems.  In  general  it  deals  with  tropospheric 
turbulence  in  terms  of  its  effects  on  the  accuracy  of  line-of- 
sight  microwave  ranging  and  angle  measuring  systems. 
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The  experimental  work  in  this  category  which  we  at  NBS  hope 
to  undertake  consists  of  three  tasks;  1)  the  study  of  "bias" 
errors  in  range  and  range  difference  (angle)  measurements, 

2)  the  systematic  study  of  atmospheric  effects  on  baseline 
system  accuracy  as  a  function  of  baseline  length,  3)  the  exten¬ 
sion  of  these  studies  from  the  present  X-band  region  (9.4  Gc/s) 
to  K-band  (22-24  Gc/s) . 

1.  Although  we  have  developed  precision  equipment  to 
measure  variations  in  the  phase-of-arrival  of  X-band  signals 
on  the  order  of  0.1  degree  (corresponding  to  a  change  in  ap¬ 
parent  range  of  .01  mm)  relative  to  an  arbitrary  reference 
value,  we  cannot  at  present  make  absolute  range  or  range  dif¬ 
ference  measurements.  Hence,  although  we  have  gathered  a  con¬ 
siderable  amount  of  data  at  the  upper  end  of  the  power  spectrum 
(from  about  1  cycle  per  hour  to  10  c/s)  we  have  not  been  able 
to  look  at  the  low  frequency  end  of  the  spectrimi,  i.e.,  to 
measure  the  "biases"  caused  by  relatively  persistent  atmos¬ 
pheric  characteristics  such  as  the  vertical  gradient  of  refrac¬ 
tive  index.  The  existing  equipment  should  be  modified  to  re¬ 
move  phase  ambiguities  so  that  the  bias  errors  and  short-term 
errors  would  be  recorded  simultaneously. 

2.  In  designing  a  baseline  tracking  system,  the  baseline 
length  is  a  parameter  of  utmost  importance.  However,  except 
for  the  experiment  performed  by  NBS  in  Hawaii  in  1956,  little 
has  been  done  to  study  the  atmospheric  errors  in  range  difference 
measurements  as  a  function  of  baseline  length.  The  Hawaii  data 
lacked  the  long-term  continuity  necessary  to  study  adequately 
the  crucial  low  frequency  end  of  the  power  spectrum  of  range 
difference  variations.  Subsequent  improvements  in  the  equip¬ 
ment  make  it  possible  to  record  range  difference  variations 
simultaneously  on  several  baselines  and  continuously  for  periods 
of  several  days.  This  should  be  done  systematically  for  base¬ 
lines  ranging  in  length  from,  say,  two  feet  (of  interest  in 
geodetic  angle-measuring  applications)  to  as  long  as  the  experi¬ 
mental  setup  (terrain,  target  distance,  etc.)  will  permit.  The 
objective  of  this  study  would  be  a  determination  of  the  expected 
rms  angular  error  versus  baseline  length  with  the  possibility 

of  determining  optimum  baseline  lengths  for  various  averaging 
times .  It  should  be  noted  that  refractometer  measurements  made 
at  each  end  of  the  baseline  may  be  used  to  reduce  substantially 
the  rms  angular  errors  on  the  longer  baselines.  Measurement  of 
the  extent  of  this  possible  reduction  is  another  objective  of 


the  program.  This  study  should  include  vertical  (tower  mounted) 
baselines  as  well  as  horizontal  baselines  to  determine  the  rela¬ 
tive  merits  of  the  two  configurations,  especially  for  measuring 
very  small  elevation  angles . 

3.  Over  the  past  several  years,  we  have  been  gradually 
accvimulating  the  necessary  equipment  to  make  range  and  range 
difference  measurements  at  frequencies  around  20  to  24  6c/s . 

Such  measurements  would  be  potentially  useful  in  determining 
the  extent  to  which  this  frequency  range  with  its  higher  phase 
resolution  and  greater  immunity  to  multipath  can  be  utilized 
in  tracking  systems.  Also,  it  is  contemplated  that  measurements 
will  be  made  on  a  frequency  in  the  range  10  to  12  Gc/s  and 
simultaneously  on  its  harmonic  which  will  lie  in  the  range  20 
to  24  Gc/s.  In  this  way  it  is  hoped  that  some  experimental 
data  on  the  dispersive  effect  of  water  vapor  can  be  obtained. 

It  is  also  hoped  that  a  refractometer  can  be  developed  for  use 
in  this  higher  frequency  range.  Such  a  refractometer  would 
have  a  smaller  cavity  and  thus  be  potentially  useful  for  smaller- 
scale  tvirbulence  studies . 
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